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Extrinsic  silicon  infrared  detectors. 
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"'Operation  and  characteristics  of  the  impurity  doped  infrared 
sensing  MOSFET  (IRFET)  are  described.  Responsi vities  of  over 
3 

10  amps/joule  have  been  observed  with  only  small  channel  width 
to  length  ratios  and  large  area  devices.  Much  higher  respon- 
sivities  would  be  possible  with  smaller  area  devices. 
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^In  the  case  of  indium-  and  gallium-doped  devices  it  has  been 
shown  that  quantum  efficiencies  in  the  range  1.0  to  10.0  % 
can  be  acheived  by  using  substrates  with  low  boron  dopings. 

The  infrared  sensing  MOSFET(IRFET)  is  a  detector,  integ¬ 
rating  element,  and  amplifier  all  combined  in  one  device  struct 
ure.  As  such  it  is  distinctly  different  than  other  types  of 
infrared  detectors  and  possesses  many  very  unique  characterist¬ 
ics.  When  compared  to  CCD  scanned  photoconductors  two  of 
these  are: 

(i)  the  responsivity  of  the  IRFET  and  uniformity  of 
responsivity  do  not  depend  upon  carrier  life¬ 
time  and  residual  impurity  concentrations,'^’1' 

(ii)  the  IRFET  has  a  nondestructive  D.C.  or  static 
memory  type  readout  as  opposed  to  the  A.C.  or 
dynamic  memory  type  of  CCD's. 

This  report  describes  in  detail  the  experimental  results 
obtained  on  gold-,  indium-,  and  gallium-doped  devices  for  use 
in  the  near,  middle,  and  far  infrared  wave  length  regions 
respectively.  Equations  describing  the  operation  of  the 
devices  have  been  compared  to  the  observed  thermal  and  optical 
response  characteristics.*-— A  good  correspondence  has  been 
found  in  all  cases.  This  allows  an  Understanding  of  the  design 
considerations  involved  in  the  MOSFET  detector  element  itself. 

An  investigation  of  the  noise  characteristics  of  the 
detector  have  been  presented  and  a  description  given  of  the 
considerations  involved  in  the  design  and  application  of  large 
scale  arrays. 

It  is  proposed  that  the  IRFET  might  be  particularly 
usefull  in  some  types  of  infrared  imaging  applications. 

Apart  from  application  of  the  MOSFET  device  structure 
as  an  infrared  detector  the  experimental  results  have  also 
demonstrated  the  usefull ness  of  the  device  structure  in 
characterizing  impurity  levels  in  silicon.  The  measurement 
of  the  thermal  and  optical  emission  characteristics  of  indium 
and  gallium  in  silicon  constitutes  to  the  author's  best 
knowledge  the  first  direct  observation  of  these  quantities. 
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TECHNICAL  REPORT  SUMMARY 


The  infrared  sensing  MOSFET  operates  on  the  principle  of  an  indirect  obser¬ 
vation  of  impurity  charge  states  in  the  surface  depletion  region  of  the  MOSFET 
device  structure.  Modulation  of  the  charge  state  of  the  impurity  center,  either 
gold,  indium  or  gallium,  in  the  space  charge  region  causes  a  change  in  the  MOSFET 
threshold  voltage  and  modulation  of  its  conductivity.  The  device  is  an  integra¬ 
ting  detector  and  must  be  periodically  reset  by  turning  the  MOSFET  off  by  accumu¬ 
lating  the  surface.  The  device  works  much  like  a  static  or  D.C.  ready  only  memory 
element.  Low  temperature  operation,  100°K,  50°K,  or  20°K  respectively  for  each 
impurity  is  required  to  avoid  thermal  emission  and  ionization  of  the  impurity 
centers.  At  low  temperatures  in  a  depletion  region  then,  only  optical  emission  of 
photoinization  is  possible  as  a  means  for  the  impurity  centers  to  change  charge 
state.  This  photoionization  is  then  the  basis  for  device  operation. 

In  operation  as  an  infrared  detector,  the  device  is  first  reset  by  accumu¬ 
lating  the  surface  and  filling  all  the  impurity  centers  with  holes  and  leaving 
them  in  the  neutral  charge  state  in  the  case  of  simple  accept  for  centers  as 
indium  and  gallium.  The  MOSFET  is  then  operated  in  the  on  conduction  state  by 
applying  an  inversion  voltage  to  the  gate.  This  inversion  voltage  and  any  back- 
gate  or  substrate  bias  forms  a  surface  depletion  region.  If  the  neutral  acceptor 
centers  subsequently  change  charge  states  from  neutral  or  negative  tne  number  of 
electrons  or  amount  of  negative  charge  in  the  inversion  layer  or  channel  must 
decrease,  resulting  in  a  lo wer  conductivity. 

The  results  presented  describe  the  fabrication,  operation  and  characteriza¬ 
tion  of  gold-,  indium-,  and  gallium-doped  devices.  Not  only  was  operation  of 
the  devices  demonstrated,  but  is  was  also  shown  that  the  MOSFET  device  structure 
possess  some  unique  abilities  in  the  characterization  of  impurity  centers  in 


silicon.  Specifically,  measurements  have  been  made  of  the  thermal  emission 

rates,  thermal  ionization  energies,  and  field-enhanced  thermal  emission  by  the 

Poole-Frenkel  effect  for  the  gallium  center. 

The  photoionization  cross  sections  of  gold,  indium,  and  gallium  in  silicon, 

which  determine  the  spectral  response  characteristics  of  the  detector  have  been 

measured.  A  strong  field  enhancement  has  been  observed  in  both  the  thermal  and 

optical  emission  from  the  gallium  center. 

It  has  been  found  that  optimum  device  performance  can  be  achieved  with  sub- 

13  14  3 

strates  with  low  boron  concentrations  of  the  order  1x10  to  1x10  /cm  and 

16  3 

indium  or  gallium  concentrations  of  the  order  5x10  /cm  .  Such  devices  will  give 

quantum  efficiencies  of  the  order  1.0  to  10.0%  and  responsivities  in  the  range 

4  9 

1x10  Amps/ Joule  for  large  area  devices  and  1x10  amps/ Joule  for  small  area 

devices. 

A  comparison  has  been  made  of  the  gallium  doped  M0SFET  characteristics  to 
results  published  for  CCD  scanned  gallium  doped  photoconductors.  While  the  quantum 
efficiency  of  the  IRFET  might  be  somewhat  lower,  on  the  other  hand,  it  has  a  much 
higher  responsivity  and  the  theoretical  signal  to  noise  ratio  under  high  flux 
conditions  is  at  least  comparable  to,  if  not  larger  than,  that  for  CCD  scanned 
photoconductors.  The  optical  responsivity  and  unformity  of  responsivity  of  the 
IRFET  is  distinctly  different  than  the  photoconductive  detectors,  in  that  the 
IRFET  characteristics  do  not  depend  upon  carrier  lifetimes  and  concentrations 
of  residual  impurities  in  the  wafers. 

As  a  consequence  of  the  unique  characteristics,  the  infrared  sensing  M0SFET 
(IRFET)  might  be  particularly  useful  in  some  large  scale  integrated  infrared 
imaging  array  applications. 
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INTRODUCTION 


DEVICE  DESIGN* 

A  design  is  described  for  an  infrared  sensing  MOSFET  (IRFET).  This  type 
of  device  has  been  designed  for  infrared  imaging  employing  a  read  only  memory 
(ROM)  array  where  individual  sense  MOSFET1 s  are  employed  at  the  intersections 
of  the  X-Y  address  lines.  The  X-Y  selection  is  provided  by  decode  circuits 
at  the  edge  of  the  array  [1].  The  resistance  of  each  IRFET  in  the  array  is 
measured  by  circuitry  external  to  the  array,  and  this  information  is  then 
converted  to  a  digital  signal  for  processing,  storage,  or  display. 

The  IRFET  is  operated  at  low  temperature  and  impurity  photoionization  [2], 
employed  to  modulate  the  conductance  of  each  IRFET.  An  individual  IRFET  is 
shown  in  Fig.  1(b).  In  the  example  illustrated,  n-channel  devices  are  employed 
and  the  substrate  is  p-type.  The  p-type  substrate  is  doped  with  both  the 
shallow  level  acceptor,  boron,  and  in  this  example,  also  by  gallium,  which  has 
an  acceptor  level  about  0.072  eV  above  the  valence  band.  At  the  temperature 
of  operation,  T  =  20  K,  the  relative  dopings  have  been  chosen  such  that  the 
gallium  acceptor  centers  are  not  ionized  [3]. 

If  the  MOSFET  has  been  originally  preset  by  turning  it  off  and  accumulating 
the  surface,  then  when  the  gate  voltage  VGG  is  applied,  the  gallium  impurity 
centers  in  the  surface  space  charge  or  depletion  region  will  not  be  ionized  at 
low  temperature  [3].  However,  the  conductance  of  the  MOSFET  is  a  function  of 
the  number  of  ionized  impurity  centers  in  the  space-charge  region  since  the 

*Based  on  material  which  has  been  published  in  the  IEEE  Trans,  on  Electron 
Devices,  Aug.  1974. 
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Fig.  5..  (a)  Capture  and  emission  at  an  impurity  center,  (b)  Infrared 


sensing  MOSFET. 
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threshold  voltage  is  a  function  of  this  space  charge  [4], [5].  If  infrared 
radiation  8  <  X  <  14  p,  illuminates  the  array,  then  the  gallium  centers  in 
the  space-charge  region  can  be  photoionized  by  the  emission  of  a  hole  to  the 
valence  band  [3]. 

Photoionization  of  the  gallium  impurity  centers  in  the  space-charge 
region  will  change  the  space  charge  and  modulate  the  resistance  of  the  IRFET. 
The  advantage  of  this  technique  for  infrared  imaging  is  the  large-scale 
integrated  ROM  arrays  can  be  constructed  with  silicon  technology.  This  is 
an  integrating  type  imaging  array;  each  device  integrates  a  fraction  of  the 
total  number  of  photons  incident  before  the  resistance  of  the  IRFET  is  sampled 
by  the  location  being  addressed. 

The  basic  physical  mechanism  employed  by  the  IRFET  is  impurity  photo¬ 
ionization.  Fig.  1(a)  shows  the  important  processes  for  an  impurity  center 
with  an  energy  level  near  the  valence  band  edge.  Three  processes  are  demon¬ 
strated:  ep^  is  the  optical  emission;  ep  is  the  thermal  emission;  and  cp  is 
the  hold  capture  [6].  Capture,  cp,  is  important  only  during  the  reset  oper¬ 
ation  of  the  device.  At  all  other  times,  the  impurity  centers  of  interest  are 
in  a  depletion  region  and  there  are  no  free  holes  present  to  be  captured  [3]. 

The  optical  emission,  ep^,  and  the  thermal  emission,  ep,  are  competitive 
processes;  the  thermal  emission  rate,  ep,  is,  however,  an  extremely  strong 
function  of  temperature  [7].  The  energy  level  of  some  impurity  centers  of 
interest  in  silicon  are  shown  in  Fig.  2.  In  this  application,  it  would 
probably  be  most  convenient  to  use  p-type  substrates  because  of  the  convenient 
placement  of  impurity  energy  levels  with  respect  to  the  valence  band.  For 
the  gold  donor  level,  the  thermal  emission  rate  is  very  small  or  the  level  is 
stable  at  the  temperature  T  =  110  K  [7],  and  for  indium  and  gallium,  the  levels 
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should  be  stable  at  50  and  20  K,  respectively.  In  the  IRFET,  the  temperature 
must  be  low  enough  to  effectively  eliminate  thermal  emission.  The  only  means 
available  for  the  impurity  center  to  change  charge  state  is  then  photoionization. 

Fig.  2  shows  the  photoionization  cross  sections,  a°Ciicj),  of  the  three 
selected  impurity  centers  in  silicon.  The  data  for  gallium  and  indium  have 
been  taken  from  Bebb  and  Chapman  [2],  while  the  data  for  the  gold  donor  level 
is  taken  from  results  obtained  by  the  author  [8].  Each  of  these  impurities 
are  most  useful  over  different  portions  of  the  infrared  spectrum. 

DEVICE  EQUATIONS 

The  basic  equations  describing  MOSFET  operation  in  the  saturation  region 

are: 

IDS  ■  (liC0)(W/L)(VGS  -  Vt)2/2 
and 

VT  =  A  +  B(Nb)1/2 

where  1^  is  the  drain  to  source  current,  y  the  effective  surface  mobility, 

CQ  the  oxide  capacitance,  W/L  the  channel  width  to  length  ratio,  VQS  the  gate 
to  source  voltage,  Vy  the  threshold  voltage,  and  A  and  B  are  parameters  which 
depend  upon  the  design  and  processing  of  the  MOSFET,  Ng  is  the  number  of 
ionized  impurity  centers  in  the  surface  space  charge  or  depletion  region. 

In  the  case  of  gold,  indium  and  gallium  doped  devices  operated  at  low 
temperatures,  Ng  can  be  time  dependent  and  is  the  combination  of  the  number 
of  shallow  level  impurity  centers,  as  for  instance  boron,  and  ionized  deep 
level  centers; 

Ng  =  +  N j Cl  -  exp(-t/x)] 


-5- 


where  all  of  the  deep  level  impurity  centers  are  neutral  or  not  ionized  at 
time  t  =  0  and  subsequently  become  ionized  by  thermal  emission  of  holes 
with  time  constant  t  =  l/e!j  ,  where  e!j  is  the  thermal  emission  rate. 

r  r 

There  are  three  distinctly  different  experimental  conditions  which  may 
exist  depending  upon  the  ratio  of  concentration  of  deep  level  impurities  to 
shallow  level  impurities  and  depending  upon  the  magnitude  of  the  excess  of 
gate  voltage  above  threshold  in  comparison  to  the  change  in  threshold  voltage 
A Vy,  as  the  impurity  centers  change  charge  state. 

Case  1.  Nj  «  and  AVy  «  (V^  -  Vj) 

If  Ny  «  N^,  then; 

AVT(t)  =  VT(t)  -  Vj  ,  where  vj  =  Vy(0) 
and 

AVy(t)  =  B[Ni/(2Na1/2)][1  -  exp(-t/T)] 

where 

AVyax  =  B[Ni/(2Na1/2)]. 


If 

AVj  <<:  (VQS  "  VT^’ 


then 

IDS  =  (  C0)(W/L)(Vgs  -  v{)  [1  -  AVT(t)] 
and  if  the  change  in  drain  to  source  current  is  defined  as, 

AI0SU)  =  "  !DS 

then; 

AInc(t)  »  (pC  )(W/L)(Vr,  -  v|)[Av!Tax 


-  AV,(t)] 


AID$(t)  =  8(Vgs  -  Vj)&V^QA  exp(-t/x) 

where 

8  =  (pCo)(W/L) 

and  the  change  in  drain  current  follows  a  simple  exponential  time  dependence. 
This  corresponds  to  the  case  of  gold-doped  devices  which  will  be  discussed 
in  more  detail  later.  The  change  in  threshold  voltage  and  drain  to  source 
current  have  been  illustrated  in  Fig.  3. 

Case  2.  Nj  « 

In  the  case  of  devices  fabricated,  on  substrates  which  are  doped  primarily 
with  the  deep  level  impurity,  Nj,  and  the  shallow  level  impurity  concentration 
is  low,  «  Nj ,  two  approaches  are  possible.  The  first  approach  is  to  deal 
directly  with  the  problem  at  hand  and  obtain  an  exact,  if  somewhat  inconvenient 
expression  for  the  time  dependence,  and  the  second  to  approximate  the  time 
dependence  by  exponential  decays.  The  indium-doped  devices  employed  are  of  the 
type  Nj  «  Nft. 

If  we  let  U(t)  =  [Igs(t)]^2  then; 

U(t>  =  (B/2)1/2[Vgs  -  VT(t)] 
and 

U(t)  =  (0/2)1/2(Vgs  -  A  -  B[Na  +  Nj (1  -  exp(-t/x))]1/2) 
and  for  times  such  that 

Nj  Cl  -  exp(-t/x)]  »  Na  , 

AUmax  =  U(0)  -  U(»)  =  (lL)1/2  -  (lL)1/2 


and  if 


AU(t)  =  AUmax[l  -  (1  -  exp(-t/  ))‘^] 


-t/x  =  ln[l  -  (1  -  AU(t)/AUmax)2]. 


This  last  expression  can  be  used  to  determine  the  time  constant,  t,  associated 
with  the  change  in  impurity  charge  state  by  plotting  on  semilogarithmic  paper. 
This  technique  has  been  used  in  the  analysis  of  gallium  doped  devices. 

Another  approach  is  to  look  directly  at  the  equation 


where 


for  times  such  that 


!DS(t)  =  e[VGS  “  VT  ‘  AVT(t)]2/2 


/T(t)  =  8(Ni)1/2[1  -  exp(t/t)]1/2 


Nj Cl  -  exp(-t/x)]  »  Nft 


ID$(t)  =  (6/2)  (VQS  -  v{)2[l  -  (AV!j!ax/(VGS  -  v|) )  (1  -  exp(-t/x)  )1/2]2 


where 


AV™aX  =  B(NX)1/2 


and  this  can  again  be  broken  down  into  two  limiting  cases.  If  in  the  special 


AVTaX  =  (VQS  -  Vj)  ; 


or  in  other  words 


!DS  =  0  ’ 


V.Y %■  >Vw; >-*  VvY^S  -•->  •  V  vVV.-VvV  W V 


«'•  *  *  .N 
\‘a  ■ 


AlD$(t)  =  (6/2) (VGS  -  vj)2(l  -  [1  -  exp(-t/T)]1/2)2  . 

This  last  expression  can  be  approximated  by  an  exponential  decay  with  a  time 
constant  of  about  one  quarter  of  the  time  constant  of  the  change  in  impurity 
change  state,  t.  If  on  the  other  hand 

AVj  «  (VgS  -  Vy)  , 

but  still 

NI  >>  NA 

then 

A!DS(t)  =  IDS(t)  -  IDS 

AIds.  (t)  =  (6/2)  (Vgs  -  Vy  )2  (ZAVyjnax  /(V^-  v/))  (1- (1-exp (-t/T))1/2) 

and  now  the  decay  in  the  current  is  approximately  exponential  with  a  time 
constant  one  half  that  of  t. 

Case  3.  AVj  «  (Vg^  -  Vy);  but  the  doping  may  vary 

In  the  situations  where  AVy,  the  change  in  threshold  voltage,  is  always 
small  in  comparison  to  the  excess  of  gate  voltage  above  threshold,  then 

AIDS  =  (tjCo)(W/l-)(VGS  " 

and 

AVy(t)  =  B(Na  +  Ni)1/2(1  -  (1  -  (Nj/(Nj  +  NA))exp(-t/x))1/2) 

then,  if  Nj  «  NA,  this  just  reduces  to  Case  1,  where  the  decay  in  the  current 
is  exponential  with  time  constant  x.  If  Nj  =  NA,  it  is  also  found  that  the 
decay  in  current  is  very  closely  approximated  by  an  exponential  with  time 


constant  t,  the  same  time  constant  as  the  change  in  impurity  charge  state. 

For  the  case  Nj  >>  N^,  it  is  found  that  the  decay  in  current  has  a  time 
constant  one  half  that  of  t,  but  just  again  like  Case  2  where  Nj  »  and 

iVT  «  <VGS  -  VT>- 

GOLD,  INDIUM,  AND  GALLIUM  DOPED  DEVICES 

Because  of  the  prior  experience  of  the  principal  investigator  with 
gold-doping  of  silicon  [8] ,  the  initial  work  was  conducted  on  gold-doping 
devices.  The  objective  of  this  work  was  to  demonstrate  that  the  infrared 
sensing  MOSFET  would  work  in  the  manner  described  in  the  original  proposal 
In  doing  so,  however,  care  was  exercised  not  to  overcompensate  the  substrate 
with  gold  which  may  cause  the  substrate  to  change  conductivity  type,  or  at 
least  become  nearly  intrinsic  and  highly  resistive  even  at  room  temperature. 

As  a  consequence  the  gold-doped  devices  are  of  the  type  Nj«N^,  or  have  an 
impurity  or  gold  concentration,  Njless  than  the  boron  concentration,  NA. 

In  this  case  the  simple  model  discussed  in  previous  section  applies. 

Subsequently  it  was  realised  that  much  higher  quantum  efficiencies  could 
be  achived  in  the  case  of  indium  and  gallium  doped  devices  by  using  substrates 
with  low  boron  concentrations. This ,  however,  result  in  a  difference  in  the 
time  constant  between  the  decay  of  the  impurity  charge  state  and  the  time 
constant  in  the  decay  of  drain  current.  In  the  case  of  indium  doped  devices  a 
numerical  solution  was  used  to  determine  this  difference  as  exactly  as  possible. 

In  the  case  of  the  gallium  doped  devices  the  approximate  solutions 

1/2 

presented  in  the  previous  section  where  employed  and  specifically  I^g  was 
plotted  for  devices  where  Nj>>  to  determine  the  time  constant  of  the 
change  in  gallium  charge  state. 

As  a  result,  while  it  might  appear  that  some  of  the  material  to  be 
presented  in  the  next  three  sections  on  each  of  the  impurities  is  repetitious 


on  the  other  hand  in  reality  each  is  being  treated  starting  for  a  different 
view  point.  This  treatment  reflects  the  historical  development  of  the  work 
and  reflects  the  realization  about  one  half  way  through  the  project  that 
many  objections  about  low  quantum  efficiency  might  be  overcome  by  using 
devices  with  The  price  in  doing  so  however  is  a  considerable  comp- 

-lication  in  the  device  model  and  equations  describing  its  operation. 

In  addition,  the  gold  doping  is  considerably  different  than  the  indium 
and  gallium  doping.  Indium  and  gallium  are  simple  acceptors,  no  amount 
of  indium  and  gallium  doping  will  cause  a  boron  doped  wafer  to  become  highly 
resistive  at  room  temperature .  Gold  on  the  other  hand  is  a  double  level 
impurity  with  both  a  deep  acceptor  and  deep  donor  level.  As  a  restilt  gold 
doping  will  cause  either  p-type  or  n-type  wafers  to  become  highly  resistive 
even  at  room  temperature .  As  a  result  gold-doped  devices  must  be  of  the  type 
NI--NA  while  indium  and  gallium  doped  devices  should  be  of  the  type  Nj»N^* 
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GOLD-DOPED  DEVICES^ 


Gold  has  been  used  most  commonly  as  a  deep  impurity  in  silicon  diodes  and 
transistors  to  increase  their  switching  speed.  An  accurate  description  of  the 
behavior  of  gold  in  silicon  is  useful  in  designing  these  devices  and  predicting 
their  characteristics;  therefore,  exhaustive  work  has  been  done  by  several 
investigators  in  order  to  describe  the  thermal  and  optical  emission  and  capture 
of  holes  and  electrons  at  gold  centers  in  bulk  silicon  [1,2,3].  This  paper 
discusses  an  extension  of  the  work  carried  out  so  far  on  the  influence  of  gold 
in  the  surface  space-charge  region  of  MOS  devices.  Our  work  is  the  first 
known  attempt  to  fully  examine  the  effects  of  infrared  radiation  on  the  gold 
impurity  centers  in  a  surface  depletion  region  behind  the  inverted  channel  of 
metal-oxide-silicon  field-effect  transistor  (MOSFET). 

As  a  deep  two-level  impurity  in  bulk  silicon,  gold  introduces  an  acceptor 
level  and  a  donor  level  in  the  bandgap  [4,5].  The  gold  center  can  be  ionized 
by  the  capture  of  holes  from,  and  the  emission  of  holes  to,  the  valence  band 
and/or  the  capture  of  electrons  from,  and  the  emission  of  electrons  to,  the 
conduction  band  (see  Fig.  1[6]).  This  ionization  will  occur  as  a  result  of 
any  sufficient  thermal  or  optical  excitation.  Within  a  depletion  region, 
however,  only  the  emission  of  carriers  from  the  centers  is  probable,  since 
there  are  no  free  holes  or  electrons  available  for  recapture  until  the  reverse- 
bias,  or  inversion  condition,  is  removed  by  accumulation.  In  addition,  the 
thermal  emission  rate  is  a  very  strong  function  of  temperature[7] ;  thus, 
thermal  emission  can  be  effectively  eliminated  by  adequate  cooling. 

Gold-doped  MOSFETs  have  been  designed  and  fabricated  in  the  University 

of  Arkansas  Solid-State  Devices  Laboratory  for  use  as  infrared  detectors 

*Th1s  section  is  based  on  material  which  has  been  published  in  the  IEEE 
Trans,  on  Electron  Devices,  Oct.  1975. 
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GOLD  IMPURITY  CENTER  IN  SILICON 


Fig.l.  The  important  carrier  emission  processes  from  the  energy  levels  of 
gold  in  silicon.  In  each  case  the  superscripts  "t"  and  "o"  denote  thermal 
and  optical  phenomena,  respectively.  ep-i  the  emission  of  holes  from  the 


gold  donor  level  to  the  valence  band,  e  is  the  emission  of  holes  from  the 

po 


gold  accepter  level  to  the  valence  band,  and  e^is  the  emission  of  electrons 
from  the  acceptor  level  to  the  conduction  band. 


(IRFETs)  [8].  The  basic  research  objective  was  to  experimentally  characterize 
the  operation  of  this  new  type  of  photon  detector  [9]  when  illuminated  by 
near-infrared  wavelength  radiation,  i.e.,  1.38  to  3.54  pm.  It  was  also  desired 
to  show  that  IRFET  channel  conductance  depends  upon  the  characteristics  pre¬ 
viously  observed  for  the  gold  center  in  bulk  silicon  by  p-n  junction  capacitance 
measurements. 

THEORETICAL  OPERATION  OF  IRFETs 

The  IRFET  is  operated  at  low  temperatures,  and  impurity  photoionization  is 
employed  to  modulate  the  channel  conductance.  In  the  schematic  diagram  of 
Fig.  2,  an  n-channel  device  is  shown.  The  p-type  substrate  is  doped  with  both 
the  shallow-level  acceptor  boron,  and  gold.  At  the  temperature  of  operation, 
the  relative  dopings  have  been  chosen  such  that  all  the  gold  centers  are  ionized. 

If  the  MOSFET  is  originally  preset  by  turning  it  off  and  accumulating  the 
surface,  then  the  gold  impurity  centers  in  the  surface  space-charge,  or  deple¬ 
tion,  region  are  ionized  by  the  capture  of  holes  (positive  charge  state)  and 
remain  filled  when  positive  gate  voltage  is  applied.  In  theory,  the  conductance 
of  the  MOSFET  channel  is  a  function  of  the  number  of  ionized  impurity  centers  in 
the  space-charge  region  since  the  threshold  voltage  is  a  function  of  this  space 
charge.  If  infrared  radiation,  e.g.,  1.77  to  3.54  ym,  illuminates  the  device, 
then  the  gold  centers  in  the  depletion  region  can  be  photoionized  or  discharged 
by  the  emission  of  holes  to  the  valence  band.  By  restricting  the  discussion 
to  the  donor-level  response,  analysis  of  the  IRFET  is  greatly  simplified.  At 
higher  photon  energies,  the  response  involves  the  combined  emission  rates  of 
the  acceptor  and  donor  levels. 

Photoionization  of  the  gold  impurity  centers  in  the  surface  depletion  region 
will  change  the  net  space  charge  and  modulate  the  conductance  of  the  IRFET. 
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GOLD  -DOPED  MOSFET  (IRFET) 

Fig. 2.  A  schematic  diagram  of  the  gold-doped,  infrared -sensing  MOSFET  (IRFET). 
Qp..  denotes  the  charge  contributed  by  electrons  in  the  surface  channel. 


This  basic  operating  principle  may  be  better  understood  by  writing  the  charge 
equation  as 


QG  +  Qss  =  '  (QCH  +  QAu  +  V  ’  (1) 

where  Qg  is  the  charge  at  the  gate  electrode,  Qss  is  the  fixed  positive  surface- 
state  charge  (which  is  dependent  on  crystalline  orientation  and  process)  at 
the  Si -Si 02  interface,  QCH  is  the  charge  contributed  by  electrons  in  the  channel, 
Q^u  is  the  charge  contributed  by  ionized  gold  in  the  oxide  or  near  the  silicon 
surface,  and  Qg  is  the  net  charge  in  the  surface  depletion  region  generated  by 
ionized  doping  impurity  atoms.  Since  all  terms  in  Eq.  (1)  except  QgH  and  Qg 
are  constant,  the  n-channel  gold-doped  MOSFET  must  obviously  lose  electrons 
from  its  inversion  layer  as  the  net  surface  space  charge  becomes  more  negative 
by  photoionization  of  the  gold  centers  due  to  infrared  illumination,  and  channel 
conductance  thus  decreases  in  order  to  maintain  charge  neutrality. 

The  basic  physical  mechanism  of  impurity  photoionization  is  also  shown  in 
Fig.  2.  The  important  microscopic  processes  for  an  impurity  center  with  an 
energy  level  near  the  valence  band  edge  of  silicon  are  optical  hold  emission 
e°,  thermal  hole  emission  e^  (which  is  negligible  at  low  temperatures),  and 
hold  capture  c^.  Capture  is  important  only  during  the  reset  operation  of  the 
IRFET  since  at  all  other  times  the  gold  impurity  centers  are  in  a  depletion 
region  where  there  are  no  free  holes  to  be  captured. 

DEVELOPMENT  OF  A  DEVICE  MODEL 

A.  Threshold  Voltage 

The  turn-on  voltage  Vj  of  a  MOSFET  is 

VT  -  VFB  ♦  24y  -  Qs/C0  ,  (2) 

where  Vpg  is  the  gate  voltage  required  to  establish  the  flat-band  condition 


(no  charge  is  induced  in  the  silicon  substrate)  and  2<j>^  -  Qg/CQ  is  the  gate 
voltage  required  to  just  bend  the  bands  through  an  amount  equal  to  twice  the 
Fermi  potential  <J>^.  The  gate  oxide  capacitance  per  unit  area,  CQ,  is  defined 
as 


Co  ”  Koxeo^ox 


where  <  is  the  relative  dielectric  constant  of  silicon  dioxide,  e  is  the 

U  A  U 

electrical  permittivity  of  free  space,  and  xn¥  is  the  thickness  of  the  oxide. 

Eq.  (2)  is  well-known  and  a  further  definition  of  terms,  or  a  more  detailed 
derivation,  may  be  found  in  several  reference  works  [10,11]. 

The  channel  conductance  of  a  MOSFET  can  be  modified  by  application  of  a 
reverse-bias  voltage  Vgg  to  the  substrate.  This  back-gate  bias  serves  to 
ionize  additional  dopant  impurities  in  the  space-charge  region  beneath  the 
inversion  layer.  If  boron  is  the  only  impurity  doping  in  the  silicon  substrate, 
back-gate  bias  affects  the  threshold  voltage  of  Eq.  (2)  through  the  term  Qg. 

The  expression  for  the  threshold  voltage  of  an  n-channel  MOS  device  thus  becomes 


VT  =  VFB 


+  2<Jy  + 


where  k  is  the  relative  dielectric  constant  of  silicon,  q  is  the  electronic 
s 

charge,  and  is  the  number  of  boron  acceptor  atoms  per  unit  volume  of  the 
substrate.  As  back-gate  bias  increases  in  absolute  value,  the  threshold 
voltage  for  a  depletion-mode,  or  normally  "ON",  device  decreases  in  absolute 
value. 


B.  Change  in  Threshold  Voltage 


the  gold  impurity  centers  (positive).  Eq.  (4)  may  be  rewritten  as 

VT  =  A  +  B  ,  (5) 

where 

A  =  VFB  +  2<j)f ,  (5a) 

8  =  ^Kseoq(^~  -"Vb"gT  /Co  (5b) 

and  Nj  is  the  effective  impurity  doping  concentration.  In  Eq.  (5),  the  Fermi 
potential  is  assumed  to  remain  approximately  constant  for  a  device  with 
substrate  bias,  independent  of  temperature. 

If  the  boron  doping  concentration  is  much  greater  than  the  gold  doping 
concentration,  and  the  cooled  (T  <  100°K)  IRFET  is  preset  in  the  dark,  then 
all  the  gold  centers  will  be  filled  with  holes,  or  positively  charged.  The 
effective  doping  concentration  is  Nj  =  NA  -  Pp,  where  Pp  is  the  total  number 
of  ionized  gold  donor  atoms.  The  gold  donor  centers  in  the  surface  depletion 
region  behind  an  inverted  channel  remain  in  their  positive  charge  state  until 
the  device  is  illuminated  with  near-infrared  radiation.  They  then  become 
neutral  with  an  average  time  constant  t  =  l/e°  [12],  and 

pD(t)  =  NTTexp(-t/x)  ,  (6) 

where  Nyy  is  the  number  of  gold  impurity  atoms  per  unit  volume  and  Pp(t)  is 
the  number  of  ionized  gold  donor  centers  as  a  function  of  illumination  time  [8]. 
After  surface  accumulation  at  t  =  0,  Eq.  (5)  becomes 

vj  =  A  +  B  \  -  Ntt  ,  (7) 

where  vj  is  the  initial  threshold  voltage  of  the  device  before  illumination. 

If  near-infrared  radiation  illuminates  the  IRFET  until  all  the  gold  centers 
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emit  their  holes  to  the  valence  band,  or  discharge,  then  =  0  and  Eq.  (5) 
becomes 

Vj  =  A  +  B  ,  (8) 

where  v|  is  the  .inal  threshold  voltage  of  the  device  with  all  the  gold  centers 
neutral . 

The  static  change  in  threshold  voltage,  AVj,  can  now  be  determined  by  sub¬ 
tracting  Eq.  (7)  from  Eq.  (8): 


avt  =  v|  -  v{  =  B(^  -  ATa  -  ntt)  , 


since  then 


(9) 


AVt  =  B(Nn/2v^)  . 


(9a) 


An  expression  for  the  illumination  time  dependence  of  the  threshold  voltage, 
Vj(t),  can  be  found  by  combining  Eqs.  (5),  (6),  (7),  and  (9): 

Vj(t)  -  V)  +  AVj(t)  =  V|  +  NT  [1  -  exp(-t/T)]  .  (10) 
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C.  Statement  of  Device  Equations 

The  fundamental  current-voltage  (I-V)  characteristics  of  a  MOSFET  are 

JDS  =  eVDS^VGS  '  VT  VDS/2^ 
in  the  linear  region  where  Vg^  -  Vy  >  and 

‘ds'8'vgs-vt>2/2  <12> 

in  the  saturation  region  where  VGS  *  VT  ^  VDS’  Eqs.  (11)  and  (12)  were  taken 
from  the  derivation  of  Carr  and  Mize  [10]  in  which  the  surface  carrier  mobility 
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is  assumed  to  be  constant,  and  IQS  is  the  drain-to-source  current  flowing  in 
the  channel,  VQS  is  the  voltage  applied  between  the  drain  and  source  ’  VGS  1S 
the  applied  gate  voltage,  Vy  is  the  threshold  voltage  of  the  MOSFET,  and 
3  =  pCQW/L,  where  p  is  the  surface  electron  mobility  in  an  n-channel  device 
and  W/L  is  the  geometric  width-to-length  ratio  of  the  channel. 

It  can  be  shown  that  for  a  static  change  in  the  threshold  voltage  of  a 
MOSFET,  the  resulting  change  in  the  drain-to-source  current,  AIDS,  is 

4Ids  =  -  S(iVT)VDS  (13) 


in  the  linear  region,  and 

AIDS  =  “  ^AVT^VGS  VT  aV^ 


(14) 


in  the  saturation  region.  Solving  Eqs.  (11)  and  (12)  for  3  in  the  linear  and 
saturation  regions,  respectively,  and  combining  the  results  with  Eqs.  (13)  and 
(14)  yields  two  different  equations  which  now  specify  the  static  change  in 
threshold  voltage  in  terms  of  experimentally  observable  I-V  characteristics: 


AVy  =  - 


AIDS^VGS  VT  ■  VDS^ 


(15) 


DS 


in  the  linear  region,  and  if  AID<-  «  1^,  then 

AWVGS  "  VT^ 


AVT“- 
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(16) 


DS 


in  the  saturation  region,  where  1^  is  the  current  before  illumination. 

A  description  of  the  current  transients  versus  time  readily  follows  by 
substituting  Eq.  (10)  for  Vy  in  Eqs.  (11)  and  (12).  Now 


=  eVDS^VGS  "  VT 


| 


and 


-5-  Ntt[1  -  exp(-t/r)]  -  Vn<./2) 

2v^  11  Ui 

— —  Ntt[1  -  exp(-t/t)])2 
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(18) 


in  the  linear  and  saturation  regions,  respectively,  where  IpS(t)  is  the  drain- 
to-source  current  in  the  IRFET  as  a  function  of  illumination  time.  Eqs.  (17) 
and  (18)  are  valid  only  if  all  the  gold  centers  are  in  their  positive  charge 
state  at  t  =  0. 

EXPERIMENTAL  METHOD 
A.  Technique  of  Measurements 

Drain  current,  IDS,  and  photocurrent,  AID<-,  were  measured  as  functions  of 
drain  voltage,  gate  voltage,  temperature,  photon  energy  and  illumination  time. 
These  results  were  compared  with  theoretical  calculations  in  order  to  evaluate 
the  accuracy  of  our  model.  Relative  doping  concentrations  were  determined 
from  IRFET  threshold  voltage  versus  back-gate  bias  measurements  before  and 
after  infrared  illumination. 

The  MOSFET  devices  were  mounted  in  a  temperature-controlled  sample  chamber, 
and  a  Bausch  and  Lomb  high-intensity  grating  monochromator  was  employed  as  the 
light  source.  A  germanium  or  silicon  filter  was  used  to  suppress  all  stray 
light  and  higher  orders  of  diffracted  light.  All  current  measurements  were 
performed  with  a  Hewlett-Packard  DC  Micro  Volt-Anmeter  425A  in  conjunction 
with  a  Houston  Instrument  Omnigraphic  2000  X-Y  Recorder.  All  voltages  were 
measured  with  Fluke  Digital  Multimeters  8000A.  An  iron-constantan  thermocouple 
voltage  was  monitored  with  a  Keithley  Electrometer  602  to  determine  temperature. 


The  n-channel  MOSFETs  were  manufactured  on  1-  to  2-9  cm  boron-doped  silicon 

wafers  of  (100)  orientation.  After  the  gate  oxidation,  gold  was  evaporated  on, 

and  diffused  for  7  minutes  at  1000°C  from,  the  back  face  of  each  slice.  The 

devices  were  then  pulled  straight  out  to  the  end  of  the  furnace  and  cooled  to 

15 

room  temperature.  The  resulting  gold  concentration  of  approximately  2x10  cm 
corresponds  to  the  results  of  previous  work  with  bulk  silicon  (p-n  junctions)  [1] 

O 

The  gate  oxides  of  the  MOSFETs  used  in  the  work  are  4000  A  thick  (estimated 
by  viewing  the  oxide  perpendicularly  under  white  light),  and  the  channels  have 
an  8.17  W/L  ratio  in  a  circular  geometry.  Each  device  was  placed  on  aT0-5 
header  and  mounted  in  a  reflective  cavity  to  provide  for  frontside  and  edge 
illumination. 

RESULTS  AND  DISCUSSION 
A.  Results 

A  typical  data  curve  using  monochromatic  illumination  is  shown  in  Fig.  3. 

In  this  case  the  drain-to-source  current  is  plotted  as  a  function  of  illumination 
time,  but  AIg<.  exhibits  similar  behavior  as  a  function  of  temperature.  The 
emission  rate  of  carriers  from  the  gold  centers  thus  can  be  determined  by  noting 
the  time  t  at  which  -Al^(t)  reaches  a  value  equal  to  1/e  of  its  maximum  value. 

The  thermal  time  constant  associated  with  emission  of  holes  from  the  gold 
donor  level,  x^,  is  plotted  versus  1000/T  (°K)  in  Fig.  4,  x^  =  A  exp(AE/kT), 
where  A  is  a  constant,  AE  is  the  energy  difference,  or  location,  of  the  level 
with  respect  to  the  valence  band,  k  is  Boltzmann’s  constant,  and  T  is  the 
temperature  in  degrees  Kelvin.  The  slope  of  the  line  now  determines  AE: 

AE  =  0.19842  Alog10xD/A( 1000/T)  (19) 


r*.  > 


DARK 


INVERTED 


ILLUMINATED 


ACCUMULATED 


-T° 

p-i 


VD$  =  +IOv 
VG  =  -2v 
VBg  =  “5v 


95°  K 


5  2.066julM 

=  3.0  x  I0,4cm'2  sec 


TIME  (SEC) 


100 


Fig. 3.  A  typical  example  of  the  time  dependent  drain  current,  1^,  in  an  IRFET 
due  to  infrared  illumination. 
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which  for  the  gold  donor  level  is  0.35  eV  [2].  This  result  is  also  shown  to 

closely  match  the  results  previously  obtained  for  the  gold  donor  level  in 

bulk  silicon  and  in  a  surface  depletion  region  by  measuring  capacitance 

transients  on  p-n  junctions  [1]  and  MOS  capacitors  [9],  respectively.  If  was 

therefore  concluded  that  a  concentration  of  gold  impurity  centers  (N,,  <  N^) 

was  observed  in  the  surface  depletion  region  under  the  MOSFET  channel. 

The  optical  emission  of  carriers  from  the  gold  center  is  observable  if 

the  temperature  is  low  enough  such  that  e°  »  et.  Then  e°  =  a°$,  where  a0  is 

2 

the  photoiomzation  cross  section  in  cm  and  <5  is  the  incident  photon  flux  in 
2 

number  per  cm  per  second.  However,  gold  is  a  double-level  impurity  in  silicon, 
and  optical  response  can  be  observed  for  both  the  donor  and  acceptor  levels  in 
the  near-infrared  wavelengths  from  1.13  to  3.54  ym. 

Fig.  5  shows  the  relative  photoionization  cross  sections  a°_^  and  a°Q  + 
respectively,  of  the  donor  level  and  the  acceptor  level  as  functions  of  photon 
energy,  i.e.  from  0.38  to  0.72  eV  and  0.58  to  0.9  eV.  The  acceptor  level 
response  was  observed  by  operating  the  IRFET  at  an  elevated  temperature  where 
ep_i  1S  very  large,  but  still  cool  enough  such  that  e^  and  e^  are  small. 

The  optical  time  constant  associated  with  emission  of  carriers  from  the  gold 
acceptor  level,  t^,  is  given  by 

TA''/(ep0t<,>  =  1^p0  +  onO,»  •  <2°> 

where  e°g  is  the  optical  emission  rate  of  holes  to  the  valence  band  (neutral 
charge  state)  and  e°^  is  the  optical  emission  rate  of  electrons  to  the  con¬ 
duction  band  (negative  charge  state).  Since  the  acceptor  level  is  near  the 
center  of  the  bandgap,  neither  e°g  nor  e°-|  can  be  assumed  negligible.  Extrap¬ 
olation  of  the  curves  indicates  that  the  threshold  ionization  energies  of  the 
gold  donor  and  acceptor  levels  are  0.35  eV  and  0.55  eV,  respectively. 


The  normalized  results  snown  in  Fig.  5  have  been  corrected  for  the  spectral 
characteristics  of  the  monochromator,  but  our  experiment  was  not  calibrated  in 
terms  of  absolute  photon  flux.  There  is  considerable  difficulty  in  determining 
the  absolute  magnitude  of  photoionization  cross  sections  because  reflections 
exist  within  the  MOSFET.  In  fact,  large  discrepancies  in  determining  optical 
flux  have  been  reported  by  several  different  laboratories  that  performed 
identical  measurements  [13].  The  results  of  our  work  for  gold  in  the  surface 
depletion  region  of  a  MOSFET  are  compared  to  the  normalized  results  obtained 
by  other  authors  [1,14,15]  for  gold  in  bulk  silicon.  A  most  likely  calibration 
point  of  10  cm  for  absolute  determinations  is  indicated.  The  response  of 
the  IRFET  to  infrared  radiation  thus  is  predicted  by  the  photoionization  cross 
sections  of  gold  in  silicon. 

Knowing  the  basic  MOSFET  equations  and  identifying  both  the  temperature 
and  photon  energy  dependencies  of  the  gold-center  emission  rates  enables  a 
rather  complete  mathematical  description  of  the  important  IRFET  characteristics. 
In  Fig.  6,  drain  current  1^  is  plotted  as  a  function  of  gate  voltage  Vg<-  in 
the  linear  region  both  before  and  after  infrared  illumination  for  a  source  to 
drain  voltage,  VDS  =  1,0  V*  The  theoretical  threshold  voltage  of  the  device 
may  be  found  by  extrapolating  the  linear  portion  of  each  curve  back  to  an 
intersection  with  the  gate  voltage  axis.  The  intercepts  in  Fig.  6  are  -16.4  \/ 
and  -14.6  V.  According  to  Eq.  (11),  VgS  =  +  1/^/2  when  ID<-  =  0;  thus 

vj  =  -16.4  -  1/2  =  -16.9  V,  =  -14.6  -  0.5  =  -15.1  V,  and  AVy  =  +1.8  V. 

Most  gold-doped  devices  tested  throughout  our  work  exhibited  a  negative 
threshold  voltage  shift.  This  result  contradicts  the  recent  experimental 
results  of  Sproul  and  Nassibian  [16]  who  noted  a  positive  shift  of  turn-on 
voltage  in  gold-doped,  n-channel  devices  and  proposed  that  the  gold  can  remove 
the  fast  interface  traps  of  continuous  energy  distribution  and  can  also  add 
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Fig.  6.  The  determination  of  IRFET  threshold  voltage  from  a  plot  of  drain  current 

versus  gate  voltage  in  the  linear  region  before  and  after  illumination. 
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acceptor  interfaces  states  very  close  to  the  valence  band  edge.  Collins 
et.  al.  [17]  suggested  that  the  observed  positive  shift  can  also  be  due  to 
fixed  negative  gold  ions  in  the  oxide  very  near  the  oxide-silicon  interface. 

On  the  other  hand,  Cagnina  and  Snow  [18]  demonstrated  the  presence  of 
positively  charged  ions  in  the  oxide  near  the  interface  by  diffusing  gold 
from  the  top  face  through  the  oxide.  When  diffusing  gold  from  the  silicon 
side,  Brotherton  [19]  showed,  by  neutron  activation  analysis,  even  higher 
concentrations  of  gold  in  the  oxide  near  the  interface  than  in  the  silicon. 

If  such  concentrations  were  producing  positively  charged  ions  in  the  gate 
oxide  of  the  MOSFETs,  they  would  cause  a  negative  shift  of  threshold  voltage 
as  observed  in  our  work.  However,  a  large  negative  threshold  voltage  was 
observed  for  all  our  devices  both  before  and  after  illumination.  Such  values 
are  inconsistent  with  (100)  orientation,  1-  to  2-&*cm  silicon  wafers;  therefore, 
this  behavior  is  possibly  due  to  a  high  Qss  concentration  and/or  some  sodium 
contamination  during  processing.  Our  MOSFETs  were  thus  depletion-mode  devices, 
but  this  anomalous  result  does  not  affect  the  outcome  of  this  paper. 

If  the  surface  potential  4>s  =  2<p^  is  assumed  approximately  equal  to  unity, 
then  Eq.  (5)  explicitly  states  that  threshold  voltage  Vy  is  a  function  of 
backgate  bias  voltage  Vgg  through  the  term  /I  -  Vgg  .  Differentiating  Eqs.  (7) 


and  (8)  with  respect  to  /l  -  Vgg  gives 
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respectively,  where  C  =  B//I  -  V gG  =  /2<seoq/Co.  Threshold  voltage  Vy  is 
plotted  as  a  function  of  A  -  Vgg  in  Fig.  7  for  the  cases  where  the  IRFET 
has  been  preset  and  all  the  gold  centers  are  positive  (Nj  =  -  Nyy)  before 

illumination,  and  where  the  IRFET  has  been  illuminated  and  all  the  gold 
centers  are  neutral  (Nj  =  N^).  By  squaring  the  ratio  of  the  slopes  of  the 
two  lines,  the  ratio  of  the  two  effective  impurity  doping  concentrations  is 
given.  Thus,  from  Eqs.  (21)  and  (22), 


Substituting  the  values  obtained  from  Fig.  7  into  Eq.  (23)  yields  Nyy  =  0.2  N^. 

If  =  lO^cm'^,  then  the  gold  doping  concentration  Nyy  =  2  x  lO^cm"^. 
Furthermore,  the  oxide  thickness  x„y  is  calculated  by  solving  Eq.  (22)  with 

U  A 

o 

the  aid  of  Eq.  (3),  yielding  «  3800  A.  The  anticipated  maximum  change  in 
threshold  voltage  with  Vgg  =  -6  V  is  found  now  by  solving  Eq.  (9),  yielding 
AVy  =  +1.64  V.  This  value  is  in  good  agreement  with  the  experimentally 
observed  AVy  =  +1.8  V  (cf.  Fig.  6). 

B.  Model  Verification 

The  experimental  I-V  characteristics  of  an  IRFET  are  shown  in  Fig.  8.  By 
solving  Eqs.  (11)  and  (12)  for  y^  and  yg  in  the  linear  and  saturation  regions, 

respectively,  and  evaluating  at  several  appropriate  points  on  the  curves,  an 

2 

average  surface  electron  mobility  of  290  cm  /V* sec  is  found.  The  average  y 
is  used  in  the  device  equations  to  model  the  IRFET.  This  result  indicates 
that  the  effective  surface  electron  mobility  in  gold-doped  n-channel  MOSFETs 
is  reduced  probably  due  to  additional  impurity  scattering  off  the  gold  sites  [20] 


DRAIN- TO- SOURCE  POTENTIAL  ,  VDS  (VOLTS) 


By  calculating  a  set  of  I-V  characteristics  and  comparing  them  to  experi¬ 
ment,  the  mathematical  model  developed  in  this  paper  can  be  shown  to  describe 
an  IRFET  with  an  accuracy  of  +30%  under  all  operating  conditions.  Surface 
carrier  mobility  is  not  constant,  however,  and  depends  on  the  effective  gate 
voltage  (Vg^  -  Vy),  as  illustrated  in  Fig.  9.  Here  the  photocurrent,  -AlD<., 
is  plotted  as  a  function  of  gate  voltage  in  the  linear  region.  Although 
the  observed  characteristics  closely  follow  the  functional  form  of  Eq.  (13), 
some  decrease  is  expected  near  threshold  and  slight  deviation  is  noted  at  high 
effective  gate  voltages  due  to  field-dependent  mobility. 

If  a  set  of  device  equations  which,  more  accurately  agree  with  the  experi¬ 
mental  DC  characteristics  of  Fig.  8  is  required,  then  empirical  expressions 
can  be  written  for  mobility  and  substituted  into  the  fundamental  integral  of 
derivation  [21].  Also,  the  simple  model  makes  use  of  a  constant  depletion- 
layer  approximation.  This  assumption  of  constant  Qg  with  respect  to  the 
gate-bias  potential  VgS  can  lead  to  error  in  the  analysis  of  MOSFET  devices, 
particularly  those  operated  in  the  saturation  region.  A  detailed  treatment 
of  this  second-order  effect  has  been  made  with  the  result  that  the  basic  device 
equations  are  even  more  extensively  modified  [22,23].  Fortunately,  most  IRFET 
designs  can  be  treated  analytically  with  the  model  developed  in  this  paper. 

In  Fig.  10,  the  photocurrent,  -AI ,  1S  shown  as  a  function  of  gate  voltage 
Vq<-  in  the  saturation  region.  The  curve  closely  follows  the  functional  form 
of  Eq.  (14)  with  the  result  that  the  linear  portion  of  the  curve  extrapolates 
to  Vgg  =  Vy  +  AVy/2  =  -16  V  when  Alpg  =  0.  AVy  =  +1.8  V,  which  is  again  in 
good  agreement  with  the  predicted  AVy  =  +1.64  V.  Since  the  most  significant 
errors  are  expected  in  the  saturation  region  of  IRFET  operation,  this  result 
lends  much  credibility  to  the  simple  model. 
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photocurrent,  -AIn„,  as  a  function 


The  IRFET  can  be  operated  in  the  linear  or  saturation  mode,  but  to  achieve 
high  responsivity  it  seems  most  desirable  to  operate  the  device  in  the 
saturation  region  where  photocurrent,  AI^  is  independent  of  drain  voltage  VD<. 
The  time  dependence  of  photocurrent,  AI^(t),  in  the  saturation  region  is 
obtained  from  Eqs.  (10)  and  (14): 

AI0$(t)  =  -  8(Vgs  -  V{)AVT[1  -  exp(-t/r)]  (24) 

if  AV-j.  «  Vgj.  -  Vy.  Also,  if  the  time  period  t  is  much  less  than  the  time 
constant  t,  then  Alp,,  changes  linearly  with  time,  or 

aIDS^  =  MDS  max  x  ’  (25) 

where  AIpS  max  is  the  maximum  static  photocurrent  at  t  8  ». 

Given  that  the  IRFET  integrates  the  incident  infrared  signal  for  a  time 
period  t  «  t,  there  will  be  a  total  input  energy  E^  =  $thu)A,  where  hw  is 
the  average  photon  energy  and  A(=  WL)  is  the  active  area  of  the  device.  At 
the  end  of  the  integration  period  t,  the  drain-to-source  current  will  have 
changed  an  amount  AIp<.(t).  The  output  power  PQut  =  |  AIDS( t)  |  VQS ,  and  the 
output  change  in  permanent  until  the  IRFET  is  reset.  AIpS  thus  can  be  read 
for  an  indefinitely  long  time.  This  consideration  suggests  the  following 
rather  unique  definition: 

Responsivity  =  Output  Power/Input  Energy  (26) 

An  example  is  given  in  Fig.  11  where  a  responsivity  of  4mW/yJ  has  been  easily 
achieved. 

Alp^(t)  depends  on  the  integration  time  t,  which  is  much  less  than  x, 
the  effective  gate  voltage  Vr<.  -  VT,  which  must  be  less  than  Vn<«,  and  the 
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Fig.  11.  An  example  for  calculating  the  responsivity  of  IRFET. 
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W/L  ratio,  which  remains  constant  with  decreasing  size.  Thus  Alp^(t)  will 
always  be  the  same  for  a  given  photon  flux  $,  regardless  of  the  area  of  the 
device.  Interestingly  enough,  the  output  signal  power  increases  directly -as 
the  drain  voltage  V^,  which  can  be  arbitrarily  large  but  less  than  the  reverse- 
breakdown  voltage  of  the  drain  junction,  while  the  input  signal  energy  decreases 
directly  as  the  area  (WL).  Much  higher  responsivities  are  obviously  obtainable 
for  smaller  devices  with  the  same  W/L  ratio. 

CONCLUSIONS 

An  analytical  model  has  been  developed  whereby  the  gold-doped,  infrared- 
sensing  MOSFET  (IRFET)  can  be  described  with  accuracy  sufficient  for  most 
engineering  applications.  The  experimental  results  leave  little  doubt  that 
the  energy  levels,  or  threshold  ionization  energies,  and  photoionization  cross 
sections  of  gold  in  a  surface  space-charge  region  correspond  to  those  previously 
observed  in  bulk  silicon.  Furthermore,  the  IRFET  behaves  according  to  pre¬ 
viously  published  design  criteria  [8].  The  infrared-sensing  MOSFET  has  a  low 

quantum  efficiency  n,  i.e.,  n  =  N^.^Wd  -  1 ,5(10~^)2(10^)1 .02(10~4)  = 

-5 

3.06  c  10  ,  where  Wrf  is  the  width  of  the  surface  depletion  region,  but  has  a 

very  high  gain  G  (or  ratio  of  the  number  of  electrons  in  the  external  circuit 
to  the  number  of  incident  photons,  both  in  the  same  sample  period),  e.g., 

G  =  AlQS(t)/(q*A)  =  4(10'5)/[1.602(10'19)3.33(1014)2.5(10'3)]  =  300  (cf.  Fig.  11), 
when  compared  to  a  photoconductive  detector  [8]. 

This  device  is  an  infrared  detector  that  has  some  very  unique  capabilities, 
namely,  a  static  read-only  memory,  a  negative  photoconductivity,  an  integrating 
function,  and  a  high  responsivity  that  is  limited  only  by  some  as  yet  undeter¬ 
mined  noise  consideration.  In  addition,  large-scale-integrated  detector  arrays 
(with  polysilicon  gate  construction  or  provisions  for  backside  illumination) 

should  be  possible  with  existing  silicon  technology. 

-40- 


Determining  a  suitable  noise  figure-of -merit  and  checking  device  uni¬ 
formity  in  array  configurations  comprise  the  remaining  work  to  be  done  with 
this  IRFET.  Even  though  only  the  near-infrared  characteristics  of  gold-doped 
IRFET  have  been  investigated  in  our  work  so  far,  there  appears  to  be  no  reason 
why  response  extended  through  the  far-infrared  wavelengths,  i.e.,  3  to  14  urn, 
is  not  possible  by  doping  the  device  with  indium  or  gallium  [8]. 
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INDIUM-DOPED  DEVICES* 

INTRODUCTION 

The  IRFET  [1-5]  which  operates  at  low  temperatures  employs  impurity 
photoionization  to  modulate  its  source  to  drain  conductance.  An  n-channel 
IRFET  is  shown  in  Fig.  1,  and  is  based  upon  a  MOSFET  device  structure  [6,7,8]. 
The  p-type  substrate  is  doped  with  both  indium  and  a  shallow  level  acceptor, 
boron.  The  indium  center  can  exist  in  two  distinctive  charge  states.  The 
neutral  charge  state  corresponds  to  the  condition  when  the  indium  center  is 
occupied  by  a  hole.  When  the  indium  center  is  ionized,  the  loss  of  this 
positive  hole  to  the  silicon  valence  band  leaves  the  indium  center  in  a  negative 
charge  state. 

The  IRFET  can  be  initially  preset  by  turning  it  off  and  accumulating  the 
surface.  In  an  n-channel  device,  this  would  correspond  to  the  application  of 
a  negative  voltage  to  the  gate  electrode,  Vg<-,  which  would  bend  the  bands  at 
the  silicon  surface  in  such  a  manner  that  holes,  majority  carriers,  would 
accumulate  there.  The  preset  operation  insures  that  all  of  the  indium  centers 
capture  a  hole  and  thus  exist  in  their  neutral  charge  state. 

The  IRFET  is  now  "biased  on"  by  the  formation  of  an  inversion  layer  of 
minority  carriers  between  the  source  and  drain  regions  and  a  resulting  current 
flow,  IpS,  between  these  regions.  For  the  device  under  consideration,  this 
corresponds  to  the  application  of  a  gate  voltage,  Vg^,  more  positive  than  the 
threshold  voltage,  Vj,  of  the  device  and  a  positive  drain  potential,  VDG.  As 
shown  in  Fig.  1,  the  application  of  these  voltages  produce  a  surface  depletion 
region  within  the  device.  If  Infrared  radiation  in  the  2.0  to  7.0  micrometer 
wavelength  range  now  illuminates  the  device,  sufficient  energy  will  be  available 

Based  on  material  which  has  been  accepted  for  publication  in  the  IEEE  Trans, 
on  Electron  Devices. 
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Fig.l.  Schematic  diagram  of  the  indium-doped  MOSFET(IRFET) .  denotes  the 
charge  contributed  by  electrons  in  the  surface  channel. 


to  photoionize  the  indium  center  by  the  emission  of  a  hole  to  the  valence 
band  [9,11].  This  photoionization  process  changes  the  net  charge  within  the 
surface  depletion  region  to  a  more  negative  quantity  since  the  indium  centers 
are  being  converted  from  a  neutral  to  a  negative  charge  state.  The  change 
in  charge  state  modulates  the  conductance  of  the  surface  inversion  channel 
between  source  and  drain  since  the  conductance  is  a  function  of  the  number  of 
ionized  impurity  centers,  Nj,  in  the  space  charge  region.  Therefore,  IDS  will 
be  continually  changing  as  the  photoionization  process  progresses  up  to  the 
point  where  virtually  all  of  the  indium  centers  are  eventually  ionized  [4,5]. 

Operation  of  the  indium-doped  silicon  infrared  sensing  MOSFET  (IRFET)  has 
been  previously  demonstrated  in  a  short  report  [5].  This  work  describes  a 
more  detailed  characterization  of  the  device  for  application  as  an  infrared 
detector  which  might  be  particularly  useful  in  the  3  to  5  micrometer  wavelength 
range. 

DEVICE  EQUATIONS 

The  two  basic  operating  regions  of  the  MOSFET  are  the  linear  and  saturation 
regions.  Equations  relating  drain  to  source  current  to  drain  to  source  voltage 
are  well  known  and  are  found  in  many  solid-state  device  books  [6,7,8].  In 
the  linear  region  of  operation,  -  Vy  >  VDS,  the  device  equation  is  expressed 
as 

>DS  -  *V0S<VGS  -  VT  -  V2>  <>> 


where  8  =  mCqW/L.  and  where  y  is  the  effective  surface  electron  mobility  in 
an  n-channel  device  and  W/L  is  the  geometric  width-to-length  ratio  of  the 
channel.  The  gate  oxide  capacitance  per  unit  area,  CQ,  can  be  calculated  from 


Co  =  Koxeo/xo 


(2) 


where  k  is  the  relative  dielectric  constant  of  silicon  dioxide  and  x„  is 
ox  o 

the  thickness  of  the  oxide. 

In  the  saturation  region  of  operation,  -  Vy  <  and  Vq<.  >  V-p  the 
device  equation  is  expressed  as 

*DS  =  S^VGS  “  V  /2 ‘ 

If  the  final  and  initial  charge  states  of  the  indium  impurity  center 
are  defined  as  before,  an  initial  and  final  drain  to  source  current  can  be 
expressed  as  [3] 

*DS  =  SVDS^VGS  "  VT  ‘  VDS^  ^ 

and 

!DS  =  6VDS^VGS  “  VT  "  VDS/2^ 


•  r 

in  the  linear  region  of  operation,  where  iL  and  il-  correspond  to  the  initial 


DS  DS 

and  final  values  of  the  drain  to  source  current. 


The  resulting  static  change  in  drain  to  source  current,  AI^,  may  now  be 
found  by  subtracting  Eq.  (4)  from  Eq.  (5).  Thus 


Al, 


(6) 


DS  =  *DS  “  TDS  =  "6(AVVDS 

where  AV-j.  is  the  maximum  change  in  threshold  voltage  or  -  vj.  Following 
a  procedure  similar  to  the  above,  the  static  change  of  drain  to  source  current 
in  the  saturation  region  becomes 


AI 


DS 


=  -B(AVt)(Vgs  -  v{  -  AVt/2)  . 


(7) 


The  negative  sign  in  Eqs.  [€)  and  (7)  indicates  that  I^  decreases  in 
magnitude  during  photoionization.  This  result  is  intuitively  reasonable  since 


the  threshold  voltage  increases  as  the  indium  centers  ionize  and  the  threshold 
voltage,  Vy,  is  defined  as  the  applied  gate  voltage  at  which  drain  to  source 
conduction  is  terminated. 

If  the  IRFET  is  held  at  the  constant  operating  temperature,  the  threshold 
voltage  may  be  written  as  [3], 

Vy  =  A  +  ,  (8) 

where  A  =  VRB  +  2<j>f,  B  =  (2<sieQq(2d)f  -  Vbg))1/2/Cq,  Nj  =  Na  +  NIn(t),  NJn(t) 
is  the  number  of  ionized  or  negative  indium  centers,  and  the  shallow  level 
acceptor  concentration.  Vy  is  the  flatband  voltage  and  <j>f  the  fermi  potential. 
If  the  neutral  charge  state  is  defined  as  the  initial  condition  and  the  negative 
charge  state  as  the  final  condition,  the  following  two  expressions  may  be 
written: 

Vj  =  A  +  B/N^  (9) 

and 

vy  -  A  +  B/N A  +  NIn  ,  (10) 

where  vj  and  Vy  represent  the  initial  and  final  values  of  the  threshold  voltage 

and  N-  is  the  total  indium  concentration, 
in 

The  static  change  in  threshold  voltage,  AVy,  can  now  be  determined  by 
subtracting  Eq.  (9)  from  Eq.  (10): 

AVy  =  v{  -  v{  =  B(vNa  +  (11) 

If  the  operating  temperature  of  the  IRFET  is  chosen  such  that  e°  »  e* 
for  a  given  incident  photon  flux,  the  increase  in  threshold  voltage  given  by 
Eq.  (11)  will  be  due  only  to  the  photoionization  of  the  indium  centers  [11]. 

The  indium  centers  will  become  negative  with  respect  to  time  during  the  photo- 
ionization  process  according  to 
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(12) 


Njn ( t )  =  NIn[l  -  exp(-t/x)]  , 

where  in  this  case  t  is  l/e°.  If,  on  the  other  hand,  thermal  emission  is 
dominant,  x  is  l/e^  [11]. 

An  expression  which  describes  the  current  transients  can  be  obtained  by 
simply  replacing  V-j.  by  Eq.  (8)  in  Eqs.  (1)  and  (3).  Thus 

IDS(t)  =  3VDS(VGS  "  VT  "  B(^A  +  NInL1  "  exP("t/T^  *  ^  “  vDS/2)  (13) 
and 

l[js(t)  =  (8/2) (VG$  -  vj  -  B[v^~  +  N'InD"-  exp(-t/x) J  -  /R^)]2.  (14) 

in  the  linear  and  saturation  regions,  respectively,  IpS(t)  is  the  drain  to 
source  current  in  the  IRFET  as  a  function  of  illumination  time. 

The  type  of  decay  that  IDS  will  exhibit  can  be  more  clearly  seen  if  we 
expand  the  square  in  Eq.  (14).  Thus 

I^s(t)  =  (6/2) [F2  -  2FB/Na  +  NIn(t)  +  2FB/N^  +  B2Nft  +  B2[NA  +  NJn(t)] 

-2B2^  vftA  +  Njn ( t )  ]  ( 

where  F  =  Vg<.  -  vj  and  Njn(t)  is  given  by  Eq.  (12). 

In  experimentally  characterizing  the  device,  this  report  will  primarily 
treat  the  saturation  region  of  operation  which  corresponds  to  the  most  desirable 
mode  of  operation  due  to  the  higher  responsivity  and  lack  of  dependence  of  the 
photocurrent  on  drain  to  source  voltage. 

DEVICE  FABRICATION  AND  MEASUREMENT  TECHNIQUES 


The  n-channel  MOSFETs  were  fabricated  upon  3  to  5  ohm-cm  indium-doped 
silicon  wafers  obtained  from  General  Diode  Corporation,  Framingham,  Massachusetts. 
The  n+  source  and  drain  regions  were  formed  by  a  phosphorous  diffusion  from  a 
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20  3 

5  x  10  /cm  "spin-on"  source.  A  three  mask  process  was  employed  in  the  fabri¬ 
cation  which  resulted  in  the  gate  oxide  being  the  same  thickness  as  the  field 

o 

oxide.  The  gate  oxides  of  the  experimental  devices  were  approximately  5400  A 
thick  and  the  channels  had  a  width  to  length  ratio  of  9.32  in  a  circular 
geometry.  The  oxide  thickness  was  verified  by  measuring  the  gate  capacitance 
on  a  Boonton  Capacitance  Meter. 

The  IRFET,  which  is  mounted  on  a  T0-5  header,  was  soldered  into  a  copper 

sample  holder  in  order  to  assure  good  thermal  conductivity.  This  sample  holder, 

which  for  thermal  measurements  completely  encompassed  the  IRFET,  was  then 

attached  to  a  model  AC-2-1 10,  open  cycle,  cryo-tip  refrigerator  manufactured 

by  Air  Products  and  Chemicals,  Inc.  The  cryogenic  refrigerator  incorporates 

a  Joule-Thomson  type  refrigerator  which  uses  compressed  nitrogen  and  hydrogen 

as  the  refrigerant.  A  radiation  shield,  which  is  cooled  to  78°K,  surrounds 

the  device  and  coldfinger.  A  stainless  steel  shroud  is  then  fitted  over  this 

-5 

assembly  and  a  vacuum  of  2  x  10  torr  is  maintained  within  the  shroud  in  order 
to  thermally  isolate  the  coldfinger  from  the  outside  world.  Temperature 
control  is  achieved  by  varying  the  hydrogen  pressure  at  the  inlet  of  the  re¬ 
frigerator. 

The  operating  temperature  of  the  IRFET  was  determined  by  monitoring  the 
resistance  of  a  calibrated  Cryogenic  Platinum  Resistor  obtained  from  Weed 
Instrument  Company,  Inc. 

For  optical  measurements,  a  Carborundum  Gas  Igniter,  or  Globar,  was  chosen 
as  the  source  of  infrared  radiation.  The  "Globar  Source"  was  operated  with  a 
110  volt  input  and  reached  an  operating  temperature  of  1500°C.  The  temperature 
of  the  globar  was  determined  with  an  Ircon  Radiation  Thermometer.  The  globar 
is  an  approximate  "blackbody  source". 


Circular  variable  filters  (CVFs)  were  employed  to  produce  the  monochromatic 
infrared  photon  flux.  These  filters,  obtained  from  Opitcal  Coating  Laboratories, 
Inc.,  are  optical  interference  coatings,  vacuum  deposited  on  circular  substrates. 
The  spectral  location  may  then  be  varied  by  physical  rotation.  The  two  90° 
segments  which  were  used  in  the  experiment  covered  the  wavelength  range  from 
2.5  to  8.0  micrometers,  have  possible  half-bandwidths  less  than  1.35%  and  have 
transmission  for  wavelengths  outside  of  the  pass-band  of  less  than  or  equal  to 
.1%. 

In  the  optical  experiments,  the  IRFET  was  mounted  so  that  it  could  be 
illuminated  with  infrared  radiation  from  the  back-side.  This  was  accomplished 
by  drilling  a  1/16  inch  hole  through  the  center  of  the  T0-5  header  before  the 
IRFET  was  mounted.  A  pictoral  representation  of  the  optical  apparatus  is  shown 
in  Fig.  2.  The  CVF  was  rotated  behind  a  2  mm  slit  with  the  angle  of  rotation 
factory  calibrated  in  terms  of  center  wavelength.  IRTRAN  II,  supplied  by 
Eastman  Kodak,  served  as  the  outer  shroud  window  of  the  cryo-tip  refrigerator. 

EXPERIMENTAL  RESULTS 

The  photocurrent  in  the  saturation  region  of  operation  is  shown  plotted  as 
a  function  of  gate  voltage  in  Fig.  3.  This  plot  follows  closely  the  functional 
form  of  Eq.  (7).  The  linear  portion  of  the  curve  is  shown  extrapolated  to 
Vqs  =  Vj  +  AV^/2  =  25  volts  when  ^Iq5  =  0.  According  to  Eq.  (3),  vj  in  the 

l  1  /2 

saturation  region  can  be  found  from  a  plot  of  (I^)  versus  VqV  where  the 
VG£  axis  intercept  directly  yields  vj.  Such  a  plot  results  in  vj  *  19  volts. 
Therefore,  AVy  =  2(25-19)  =  12  volts. 

The  values  obtained  for  vj  and  v|  vary  somewhat  depending  upon  whether 
they  were  calculated  in  the  linear  or  saturation  regions.  This  is  probably 
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Fig.  2.  The  infrared  source  and  optics  employed  in  the  optical  characterization. 
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Fig. 3.  The  photocurrent,  Al^  as  a  0f  gate  voltage  in  the  saturation 
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due  to  the  fact  that  the  simple  model  developed  makes  use  of  a  constant 
depletion-layer  width  approximation. 


In  order  to  calculate  the  theoretical  change  in  threshold  voltage,  Eq.  (11) 
requires  a  knowledge  of  the  doping  concentrations  of  indium  and  the  residual 
boron  impurities.  If  the  surface  potential,  <|>  =  2<j>f,  is  assumed  to  be  approx¬ 


imately  equal  to  1  Volt,  Eq.  (8)  states  that  Vy  is  a  function  of  VgG  through 


the  term  (1  -  VgG) 


1/2 


Differentiating  Eqs.  (9)  and  (10)  with  respect  to 


(1  -  VBG)1/2  yields 


dV 


d[(i  -  vBG)-] 


I _ C(N  )1/2 

,  J/2-, 


(16) 


and 


dV. 


d[(l 


T  T72T=c(NA  +  fW1/2 


Vn ,),,c] 


(17) 


BG' 


respectively,  where  C  =  (2<Sl- eQq) 1  . 

1  /2 

In  Fig.  4  a  plot  of  Vy  versus  (1  -  VgG)  '  has  been  made  for  the  cases 
where  the  indium  centers  are  neutral  (dashed  line)  and  negative  (solid  line). 
From  an  examination  of  Eqs.  (16)  and  (17),  the  effective  dopings  for  both 
charge  states  can  be  determined  from  the  slopes  of  the  lines  in  Fig.  4. 

O 

Utilizing  an  oxide  thickness  of  5400  A  it  was  determined  that  +  NJn  = 

8.92  x  10^/cm^  and  =  4.95  x  10^/cm^.  Then,  with  an  applied  back-gate 
bias  of  -2.5  volts,  Eq.  (11)  becomes  AVy  =  1.66  x  lO"^^  +  -  *4^  = 

12.02  volts.  This  result  matches  the  experimental  value  of  AVy  determined 
in  the  saturation  region  almost  identically.  Note  in  this  case,  NJn  »  N^, 
or  the  indium  impurity  concentration  exceeds  the  shallow  level  acceptor 
concentration,  N^.  The  simplifying  assumption  NyT  «  N^  used  previously 
[1,2,3]  cannot  be  applied  in  the  analysis  of  device  operation  to  be  presented 


A  typical  plot  of  Ip,,  in  the  saturation  region  versus  time  as  a  result 
of  the  thermal  decay  of  the  indium  impurity  centers  is  shown  in  Fig.  5.  At 
a  time  before  t  =  0,  an  accumulation  voltage  was  applied  to  the  IRFET  to 
insure  that  all  of  the  indium  centers  would  be  in  their  initial  neutral  charge 
state. 

A  careful  study  of  the  resulting  curve  revealed  that  IpS(t)  obeyed  a 

nearly  perfect  exponential  decay  for  t  <  tt  where  tt  represents  the  time 

XDS  xDS 

constant  of  the  IpG(t)  decay.  Therefore,  for  the  sake  of  experimental  con¬ 
venience  IpS(t)  was  approximated  by 

In<.(t)  =  S  exp(-t/tT  )  +  T  (18) 

XDS 

where  T  =  IDS («>)  and  S  +  T  =  IpS(0).  The  time  constant  of  the  current  decay 
can  then  be  determined  by  noting  the  time  at  which  IpS(t)  has  reached  a  value 
of  exp(-l)  of  the  difference  between  its  maximum  and  IDS(°°). 

In  order  to  determine  the  time  constant  of  the  indium  decay,  Eq.  (15)  can 
be  utilized  to  solve  for  t  for  a  given  plot  of  IpS(t).  The  IRFETs  which  were 
employed  in  the  experimental  measurements,  had  the  following  representative 
parameters,  tQX  =  5400  A,  2<f>f  =  1.0,  VgG  =  -3  volts,  Vy  =  +19  volts,  VGS  =  +26 
volts,  =  4.95  x  10^4/cm3,  and  Njn(°°)  =  8.425  x  10^®/cm3,  at  an  operating 
temperature  of  52°K.  By  employing  these  parameters,  F  =  7,  B  =  1.78  x  10”^, 
and  8  =  6.82  x  10"®.  Incorporating  these  constants  in  Eq.  (15)  yields  the 
relationships 

IDS(t)  =  4.63  x  10"4  +  9.1  x  10"4[l-exp(-t/x)]  -  1.26  x  10"3[l-.945exp(-t/T)]1/2 

(19) 

Eq.  (19),  which  will  be  utilized  in  interpreting  the  experimental  data,  is 
valid  between  0  <  t/x  <_  .496  since  VGg  <  Vy  at  t/x  >_  .496.  However,  since 


an  experimental  observation  of  IDS(t)  was  found  to  be  approximated  by  Eq.  (18), 


a  relationship  is  suggested  between  tt  in  Eq.  (18)  and  t,  the  time  constant 
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of  the  indium  decay,  in  Eq.  (19).  This  relationship  may  be  expressed  as 


1/t  =  a/Tj 


DS 


(20) 


whe  re 


S  exp(-t/xT  )  +  T  =  4.63  x  10"^  +  9.1  x  10~^[1 -exp(-at/xT  )] 
XDS  DS 


1.26  x  10"3[l-.945exp(-ctt/xT  )]1/2 

aDS 


(21) 


and  where  a  is  a  constant.  The  constants  S  and  T,  which  were  defined  in  Eq.  (18), 
were  evaluated  assuming  the  same  initial  and  final  values  of  IDS(t)  as  given  by 
Eq.  (19).  A  computer  numerical  integration  routine  was  then  employed  to  find 
an  a  such  that  the  mean  square  error  between  the  two  functions  in  Eq.  (21) 
would  be  minimized  for  t/x  _<  .496,  the  limiting  value  for  which  Eq.  (19)  is 
valid.  The  resulting  a  which  makes  Eq.  (21)  true,  was  found  to  be  .11446. 
Therefore,  Eq.  (20)  becomes 


x  =  8,7t 


t 
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(22) 


Within  the  integration  limits  which  were  employed,  this  value  of  a  matched 
the  two  expressions  in  Eq.  (21)  with  an  error  of  less  than  2%. 

The  indium  decay  time  constant  may  now  be  obtained  by  determining  tt 

uS 

from  a  plot  of  I^U),  as  in  Fig.  5,  and  by  the  use  of  Eq.  (22).  This 

procedure  was  carried  out  for  different  temperatures  of  operation.  The 

resulting  time  constants  were  then  plotted  on  a  log  scale  versus  1000/T(°K) 

as  shown  in  Fig.  6.  This  difference  in  the  time  constant  of  the  current 

decay,  xT  ,  and  the  time  constant  of  the  change  in  indium  charge  state  x, 
:DS 


was  incorrectly  leftout  in  previous  reports  [4,5],  the  time  constant  shown 


there  for  indium-doped  devices  with  NT„  »  N.  is  xT  not  x. 

In  A  !DS  ! 

It  has  been  determined  that  the  time  constant  associated  with  the  ioni¬ 
zation  of  an  impurity  level  can  be  expressed  as  [9,11] 


x  =  A1T"1'5exp(AE/kT)  , 


(23) 


where  is  a  constant,  AE  represents  the  location  of  the  impurity  level  with 
respect  to  the  valence  band,  k  is  Boltzmann's  constant,  and  T  is  the  tempera¬ 
ture  in  degrees  Kelvin.  A  straight  line  drawn  through  the  data  points 
employing  a  least  mean  square  error  technique  then  allows  a  calculation  of 
the  ionization  energy,  AE,  of  indium  in  silicon. 


Aloggx  +  1 .5AlogeT 
AE  =  1 1 .6047A(  1000/T)  =  0J67eV  • 


(24) 


This  value  corresponds  to  the  value  obtained  by  other  techniques  for  the  energy  ; 

level  of  indium  in  silicon  [12-17]. 

An  operating  temperature  less  than  45°K  was  chosen  for  the  optical  measure¬ 
ments  so  that  the  thermal  emission  rate  would  be  very  small.  With  the  globar 
on,  a  -40  to  +40  volt  square  wave  was  applied  to  the  gate  electrode  which 
caused  the  silicon  surface  to  alternate  between  strong  accumulation  and  strong 
inversion,  respectively.  As  in  the  thermal  measurements,  Vg,.  was  chosen  such 
that  the  transistor  would  be  operating  in  the  saturation  region  during  inversion 
and  Ip<-  was  monitored.  A  typical  resulting  trace  is  shown  in  Fig.  7. 

The  illuminated  decay  of  1^  is  caused  by  two  separate  hole  emission 
processes;  the  optical  hole  emission  initiated  by  the  incident  1500°C  photon 
flux,  e°,  and  the  optical  hole  emission  initiated  by  the  25°C  "Background" 

flux  incident  upon  the  device,  e^.  The  IRFET  is  an  integrating  detector  and  I 
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TIME  (milliseconds) 


43.2°  K  Vnc  =  *7  volts  Vo*  =  0.0  volts 


integrates  the  background  flux  in  the  wavelength  range  out  to  8  micrometers  as 
well  as  the  monochromatic  signal  from  the  1500°C  source  passing  through  the 
filter.  The  total  hole  emission  rate,  ep,  for  the  illuminated  decay  can  now 
be  expressed  as  [11] 


e 


P 


(25) 


Recalling  that  t 


1/ep,  the  following  expression  may  be  written  for 


1/t 


1/TP  =  1/TP 


1/T 


B  • 


(26) 


With  the  use  of  Eq.  (26)  and  Fig.  7,  the  photoionization  cross  section  can  be 
determined.  The  results  shown  in  Fig.  8  have  been  corrected  for  the  spectral 
characteristics  of  the  globar,  transmission  of  the  CVF,  and  the  bandwidth  of 
the  CVF,  but  the  experiment  was  not  calibrated  in  terms  of  absolute  photon 
flux.  Considerable  difficulty  exists  in  determining  the  absolute  magnitude  of 
photon  flux  due  to  the  reflections  which  exist  within  the  IRFET  [3].  The 
photoionization  cross  section  determined  by  other  authors  [16]  from  measure¬ 
ments  on  an  indium-doped  photoconductor  is  also  shown  in  Fig.  8.  In  Fig.  8 
it  can  be  seen  that  the  general  shapes  of  the  curves  are  similar  for  both 
segments  1  and  2  of  the  CVF.  The  photoionization  cross  section  of  indium  in 
silicon  [12-16]  thus  determines  the  response  of  the  IRFET  to  incident  infrared 
radiation. 

When  Eq.  (26)  was  used  to  calculate  the  photoionization  cross  section, 

xT  ,  not  x,  was  the  parameter  employed.  However,  since  the  presented  data 
*DS 

is  in  terms  of  a  normalized  photoionization  cross  section,  the  constant  of 

proportionality  which  exists  between  xT  and  x  will  have  no  effect  upon  a 
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normalized  representation. 


NORMALIZED  PHOTOIONIZATION  CROSS-SE 


Wavelength  (microns) 


Fig.  8.  The  exper mentally  determined  photoinzation  cross  section  of  indium  in 
silicon  (solid  line)  superimposed  upon  the  result  of  Pines  and  Baron  [16] . 


In  order  to  achieve  a  high  responsivity,  the  IRFET  should  be  operated 
in  the  saturation  retion  where  the  photocurrent,  AI^,.,  is  independent  of  the 
drain  to  source  voltage.  Since  the  IRFET  integrates  the  incident  infrared 
photon  flux,  an  expression  for  the  incident  input  energy,  E^,  becomes 

Ein  =  4>thvA  ,  (27) 


where  t  is  the  length  of  the  integration  period,  A  is  the  active  area  of  the 
device,  4>  is  the  photon  flux  and  hv  is  the  photon  energy. 

At  the  end  of  an  integration  period  the  drain  to  source  current  will 
have  changed  by  an  amount  equal  to  AIDS(t).  The  fact  that  AIDS(t)  will  remain 
constant  until  reset  and  thus  can  be  "read"  for  an  indefinitely  long  period  of 
time  following  integration  allows  the  following  definition  of  IRFET  respons¬ 
ivity  to  be  proposed  [3] : 


RESPONSIVITY  = 


Photocurrent 
Input  Energy 


(28) 


An  example  of  calculating  the  responsivity  of  an  IRFET  is  shown  in  Fig.  9  [3]. 

The  photoionization  cross  section  of  indium  at  2.066  ym  was  taken  to  be 
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8.4  x  10  cm  [12,13,14].  In  doing  so  we  have  used  the  results  of  Bebb  and 
Chapman  [13]  for  the  photoionization  cross  section  rather  than  the  smaller 
value  reported  by  Messenger  and  Blakemore  [15].  The  incident  photon  flux  was 
then  determined  from  the  time  constant  of  the  optical  decay  and  photoionization 
cross  section  as  follows: 

.,nv  .  _hy_  =  h(2.998xl08)/(2.066xl0~6)  = 

T°a°  (ft. 7x3. 35) (8. 4x10” 17) 

<}>*nv  =  3.93  x  10”5  watts/cm^  ,  <J>  =  4.09  x  10^/cm^sec  (29) 


Values  of  around  86  milliwatts/microjoule  and  1.6  mil li amps/microwatt  were 

obtained  for  the  responsivity.  Much  higher  responsivities  are  expected  at 
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IRFET  RESPONSIVITY 


INDIUM  DOPED 
T  =  42°  K 

W/L  =  9.32  G  CH 

AREA  =  1.362  xIO'2  cm2 
X0  =  5400  A 


VDS? 

D  1  [  IDS  ~  !aIDS 


S  4 


A  =  2.066  microns 


“  VDS  =  +18  volts 
VBg  =  "3  volts 
VGS  =  +26vo!ts 


INPUT  ENERGY  AND  POWER 
CT°  =  8.4  x  I0“i7cm2 

Z °  s  3.35  sec,  5  s  4.09  x  10 14  /cm2  see 
j-cs 

P  =5.35  X  10"7  watts 

Ejn  (t=U.o/|0)  =  1.79  x  I0"7  Joules 
J-DS 


CHANGE  IN  OUTPUT  CURRENT  t  22  t£/IO 

■^DS 


AI 


DS 


=  .86mA 


RESPONSIVITY 


PHOTOCURRENT 
INPUT  ENERGY 


4.8  mi  Hi  amps 
microjoule 

Fig.  9.  An  example  for  calculating  the  responsivity  of  an  indium-doped  IRFET. 
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longer  wavelengths  since  the  magnitude  of  the  photoionization  cross  section 
peaks  at  approximately  4.6  micrometers.  In  addition,  since  the  photocurrent 
depends,  in  saturation,  only  upon  the  integration  time,  the  effective  gate 
voltage,  and  the  W/L  ratio,  it  will  always  be  the  same  for  a  given  incident 
photon  flux.  Therefore,  much  higher  responsivity  can  be  obtained  for  smaller 
devices  with  a  constant  W/L  ratio  since  the  input  energy  (or  input  power) 
decreases  as  the  area  (W* L)  decreases. 

APPLICATION  CONSIDERATIONS 

Operating  temperatures  of  less  than  50°K  are  required  to  insure  the 
thermal  emission  rate  of  the  indium  center  is  orders  of  magnitude  smaller 
than  the  optical  emission  rate  when  the  device  is  employed  as  an  infrared 
detector.  The  lower  limit  on  the  temperature  of  operation  is  determined  by 
the  boron  concentration  in  the  substrate  and  the  requirement  that  sufficient 
holes  are  present  in  the  substrate  to  insure  a  fast  reset  operation  where 
Treset  =  anc*»  cp’  the  caPture  coefficient  of  the  indium  center 

and,  p,  the  hold  concentration  [11,9],  For  the  devices  described  here, 

Njn  »  N^,  however,  while  deionization  of  the  indium  in  the  substrate  was 
almost  complete,  sufficient  boron  ionization  exists  such  that  the  reset  time 
in  the  20°K  to  50°K  range  was  in  the  microscecond  or  sub-microsecond  time 
scale  and  could  not  be  conveniently  determined. 

Much  higher  responsi vities  could  be  achieved  with  small  area  devices, 
on  the  other  hand,  it  has  been  shown  that  under  shot  noise  or  background 
limited  operation  the  maximum  obtainable  signal  to  noise  ratio  depends  only 
on  the  surface  charge  storage  density  or  total  number  of  indium  centers 
ionized  in  each  device,  S/N  =  (NjnWjA/2)^,  where  Wd  is  the  depletion  region 
width  and  A  the  area  [18].  Small  areas  result  in  small  S/N  ratios.  In  the 
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case  of  gold-doped  devices  operating  at  78°K  it  has  been  found  that  shot  noise 
or  background  limited  rather  than  1/f  noise  limited  operation  can  be  achieved 
for  integration  of  less  than  10  seconds  [18].  The  dependence  of  the  maximum 
obtainable  signal  to  noise  ratio  on  the  square  root  of  the  total  number  of 
stored  charges  is  a  common  characteristic  of  all  integrating  semiconductor 
infrared  detectors  and  applies  equally  as  well  to  the  CCD  scanned  photoconductor 
arrays  where  the  charge  is  integrated  and  stored  in  CCD  potential  wells  [16, 
19,20,21].  Since  the  amount  of  charge  which  can  be  stored  in  any  type  of  MOS 
surface  device  whether  MOSFET  or  CCD  is  of  the  same  order  of  magnitude,  the 
selection  of  the  particular  type  of  infrared  detector  should  be  made  on  the 
basis  of  other  considerations. 

CONCLUSIONS 

The  characteristics  of  the  indium-doped  silicon  infrared  sensing  MOSFET 
(IRFET)  operating  in  the  saturation  region  have  been  investigated  and  described 
in  some  detail.  Again  it  has  been  found  that  the  device  operation  can  be 
described  on  the  basis  of  a  simple  model  [1]  and  that  the  results  obtained 
on  gold-doped  devices  in  the  near  infrared  [3]  can  in  fact  be  extended  to 
shallower  level  impurities,  lower  temperatures,  and  longer  wavelengths.  This 
type  of  infrared  detector  might  be  particularly  useful  for  application  in 
large  scale  integrated  infrared  imaging  arrays  operating  in  the  middle 
infrared  or  3.0  to  5.0  pm  wavelength  range. 

The  IRFET  detector  has  some  unique  characteristics  when  compared  to 
alternate  techniques  involving  CCD  scanned  arrays  of  extrinsic  silicon 
photoconductive  detectors  [16,19,20,21].  In  the  case  of  the  IRFET,  the 
detector,  storage  or  memory  element,  and  the  read-out  and  output  amplifier 
have  essentially  been  integrated  into  a  single  device.  In  addition,  IRFET 
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GALLIUM-DOPED  DEVICES 


INTRODUCTION 


The  infrared  sensing  MOSFET  (IRFET)  with  gold  and  indium  doping  has  been 

demonstrated  as  an  infrared  detector  with  some  unique  characteristics  for  use 

in  the  2-5y  and  5-8y  infrared  regions  respectively  [1-7].  This  report  deals 

with  the  characteristics  of  gallium-doped  MOSFETs  for  use  in  the  far  infrared 

8  to  14  ym  wavelength  range.  Each  particular  infrared  region  has  associated 

with  it  particular  characteristics,  advantages,  and  problems.  Application  of 

gallium-doped  silicon  MOSFET  device  structures  in  the  far  infrared  is  no 

exception.  In  this  case,  the  background  in  a  300°K  blackbody  source  emitting 
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a  photon  flux  of  the  order  2  x  10  /cm  sec.  For  imaging,  to  detect  a  small 
temperature  variation  above  the  background  then  requires  stringent  control 
of  the  gallium  uniformity. 

The  operation  principle  of  the  IRFET  bases  upon  the  photoionization  pro¬ 
cess  of  a  deep  level  impurity  (gallium)  in  the  depletion  region  of  a  MOSFET 
at  a  temperature  below  the  freeze-out  temperature  of  the  gallium  center.  Upon 
irradiation,  the  Ga  centers  change  charge  state  from  neutral  to  negative  with 
a  time  constant  t.  This  changes  the  threshold  voltage  which  in  turn  changes 
the  drain  to  source  current  and  is  detected  by  an  external  circuit.  At  the 
end  of  the  exposure  period,  it  is  necessary  to  place  the  Ga  centers  in  the 
neutral  state  for  the  next  integration.  This  is  achieved  by  applying  an 
accumulation  voltage  to  the  gate  and  by  doping  the  IRFET  with  another  shallow 
level  impurity  (boron).  The  shallow  level  borons  assure  that  sufficient  holes 
are  available  to  neutralize  the  gallium  center  during  the  reset  period. 

The  MOSFET  structure  also  provides  a  convenient  tool  for  the  character¬ 
ization  of  shallow  level  impurity  centers  in  silicon  [8,9],  particularly  in 


•  v  -  *  ■’  ••  ■r.  •'  •*.  ■*.  •*. 


•• . 
-•v 


the  case  where  there  is  little  or  no  other  shallow  level  impurities  present 
and  the  substrate  resistance  becomes  very  high  at  low  temperatures. 


Equations  describing  the  operation  of  the  device  for  the  case  when  the 
density  of  deep  level  centers  is  small  compared  to  shallow  level  centers  and 
for  the  case  when  the  density  of  deep  level  centers  is  large  compared  to  the 
shallow  level  centers  are  presented.  These  situations  correspond  to  doping 
levels  used  to  characterize  the  gallium  center  and  to  doping  levels  used  to 
achieve  high  sensitivity  when  the  MOSFET  is  used  as  a  detector. 

Experimental  verification  of  the  device  operation  including  character¬ 
ization  of  the  Gallium  centers  as  functions  of  temperature,  backgate  bias, 
and  photon  energy  is  also  presented. 

DEVICE  EQUATIONS 

To  obtain  large  IQS  and  to  avoid  dependence  on  VQS,  the  device  is  oper¬ 
ated  in  the  saturation  region.  The  basic  equations  describing  MOSFET  operation 
in  the  saturation  regions  are  [10-12]: 

IQS(t)  =  (S/2)[Vgs  -  VT(t)]2  (1) 

and 

VT(t)  =  A  +  B(1  -  VBQ)1/2[NB(t)]1/2  (2) 

where  6  =  yC0W/L,  y  is  the  effective  surface  mobility,  CQ  the  oxide  capacitance 
per  unit  area,  W/L  the  channel  width  to  length  ratio,  IQS  the  drain  to  source 
current,  Vg<.  the  gate  to  source  voltage,  Vj  the  threshold  voltage,  A  and  B  are 
parameters  which  depend  upon  the  design  and  processing  of  the  MISFET,  Ng  the 
number  of  ionized  impurity  centers  in  the  depletion  region  and  VBG  is  the 
backgate  bias.  NR  can  be  time  dependent  and  is  the  combination  of  the  number 


of  ionized  shallow  level  impurity  centers  (boron)  and  ionized  deep  level 
impurity  centers  (gallium),  i.e., 

NB(t)  =  Nft  +  N r [1  -  exp(-t/x)]  (3) 

where  and  Nj  are  the  boron  and  gallium  concentration  respectively.  In 
Eq.  (3),  all  the  boron  centers  are  assumed  ionized  and  all  the  gallium  centers 
are  arranged  neutral  at  time  t  =  0  and  subsequently  become  ionized  by  thermal 
and/or  optical  hole  emission  with  a  time  constant  x  =  (et  +  e°)~^  [8],  where 

r  r 

e*  and  e°  are  the  thermal  and  optical  hole  emission  rate.  In  reality,  because 

r  r 

the  difference  in  the  boron  and  gallium  ionization  energies  is  small,  a  portion 
of  the  boron  centers  in  the  depletion  region  are  not  ionized  at  t  =  0,  and  will 
ionize  at  a  rate  1/x..  However,  for  thermal  measurements,  x  »  x^,  i.e.,  in 
the  practical  sense  all  the  shallow  impurity  centers  can  be  assumed  ionized. 

For  photoionization  measurements,  the  temperature  must  be  arranged  such  that 
for  the  gallium  centers  e^  «  e°.  Under  this  condition,  the  boron  thermal 
hold  emission  rate  can  be  comparable  to  the  gallium  optical  hold  emission  rate. 

In  other  words,  accurate  determination  of  the  Ga  photoionization  cross  section 
is  difficult  unless  the  boron  centers  are  first  ionized  with  photon  flux  having 
energy  less  than  the  Ga  threshold  energy.  For  convenience,  in  the  following 
discussion,  we  shall  assume  are  all  ionized  in  the  depletion  region  at 
t  =  0. 

There  are  two  distinctly  different  experimental  conditions  which  may 
exist  depending  upon  the  ratio  of  deep  to  shallow  level  impurities  concentration. 

A.  Case  1 :  Nj  « 

Defining  AV^.(t)  =  Vy(t)  -  V^(0),  from  Eq.  (13)  we  obtain 


¥*>  =  AVT,n,ax[1  -  exP(‘t/T)] 
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Since  NT  «  Nft,  AVT  _  will  be  small.  For  any  reasonable  Vrc,  AVT 

I  A  T,max  bS  T,max 

«  VGS  ‘  VT  (0).  Eq.  (1)  becomes 

IDS(t)  =  (B/2)[Vgs  -  VT(0)][VG$  -  VT(0)  -MT(t)]  .  (6) 

If  the  change  in  drain  to  source  current  is  defined  as 

AIDS(t>  "  IDS(t>  -  W>  <7> 


then 

Al DS ( t )  =  (6/2)(Vgs  -  VT(0)]6VT>rnaxexp(-t/T)  .  (8) 

Equation  (8)  indicates  that  Alg<-(t)  varies  exponsentially  in  time  with 
a  time  constant  r.  The  thermal  and  optical  emission  rates  of  the  deep  level 
impurity  centers  are  readily  obtained  by  measuring  the  change  in  drain  to 
source  current  of  the  MOSFET,  provided  the  shallow  and  deep  levels  are  ade¬ 
quately  separated.  Although  this  case  provides  a  convenient  means  of 
characterizing  the  centers,  as  an  infrared  detector  or  emaging  element  how¬ 
ever,  for  high  responsivity  AI^  is  to  be  as  large  as  possible.  In  view  of 
Eqs.  (5)  and  (8),  large  1^  can  be  obtained  if  Nj  is  made  much  larger  than 
This  is  discussed  in  the  following  section. 

B.  Case  2:  Nj  >> 


In  this  case,  AI^  is  not  an  exponentially  decaying  function.  From  Eqs. 
(1 )  to  (3) ,  we  obtain 

AV-j-(t)  =  B ( 1  -  Vbg)1/2{[Na  +  Nj ( 1  -  e"t/T)]1/2  -  (Na)1/2> 
and 
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AIDS(t)  =  (8/2){[Vgs  -  A  -  BO  -  Vbg)1/2[Na  +  NjO  -  e't/x)1/2]2 

-  [VGS  -  A  -  B(1  -  Vbg)1/2(Na  +  Nj)172]2} 

It  is  informative  to  obtain  AInc  which  gives  direct  merit  of  the  device 

OS, max  3 

and  is  defined  as 

AI0S,max  =  W0)  " 

=  (8/2){[Vg$  -  A  -  8(1  -  Vbg)1/2(Na)1/2]2  (11) 

-  [VGS  -  A  -  B(1  -  Vbg)1/2(Na  +  Nj)172]2}. 

From  Eq.  (11)  for  a  given  VG<-,  a  larger  AIDS  max  can  be  obtained  if  IqS(co) 
is  set  to  zero  which  occurs  at  VG<-  =  Vy(“)  =  A  +  B(1  -  vbq1//2(Na  +  N^)172. 

If  (1  -  VgG)^72(NA  +  Nj)^72  is  sufficiently  large,  cutoff  can  occur  before 
t  =  »,  i.e.,  the  IRFET  response  time  can  be  reduced  by  increasing  VBG  and/or 
Nj.  Under  the  condition  that  A  +  B(1  -  VBG)172(NA  +  N^)172  >  V^, 

MDS,™x  -  (B/2)CVgs  -  ft  -  BO  -  Vbg)'/2(Na),/2]2  .  (12) 

In  Eq.  (2),  AID<-  max  can  be  further  increased  by  keeping  NA  small  and  by 
increasing  VG<..  The  smallest  NA  is  determined  by  the  reset  time  of  the  IRFET. 

At  20°K  for  a  reset  time  of  microsecond,  NA  =10  cm  .  At  the  present  state 
of  the  art,  A  can  be  made  less  than  a  volt  [10].  For  large  VGS  and  small  NA, 
from  Eq.  (12)  we  have 

4IDS.max  ‘  •  <13> 

For  the  MOSFET  to  reamin  in  the  saturation  region,  we  must  have 
Vqs  £  +  Vj  which  serves  as  an  upperbound  for  Vgg.  Since  the  largest  V^g 

is  limited  by  the  breakdown  voltage  which  decrease  as  Nj  increases,  large  AIpg 
max  can  be  obtained  by  judiciously  by  chossing  N^and  VBG  to  satisfy  the  condition 
A+BCl-V^,)172^)172;;  VGg.  for  y*  600cm2/v.sec 


XVIII ‘so 


and 


X  =  1  KA,  W/L  =  20,  A  =  0.6V,  N.  =  1013cm"3  and  Vrc  =  10V,  AInc 
o  ’A  GS  DS.max 

AVj  maxC=  V-j-(°°)  -  VT(0)]  computed  from  Eqs.  (9)  and  (11)  are  presented  in 

Fig.  1.  Saturation  of  AIp<.  is  due  to  cutoff  of  the  MOSFET.  From  this 

figure  for  NT  =  1016cm‘3,  \lDr  =  -20 V,  a  AInc  of  16.5  mA  and  AVT  ,  of 

i  du  u^,max  I ,max 

7 V  are  obtained. 

To  illustrate  the  operation  of  the  device,  experimental  values  of  AV^ 

I  /O 

versus  (1  -  VgG)  '  and  the  change  in  drain  to  source  current  with  and 

without  radiation  for  a  device  with  somewhat  different  parameters  are  shown 

in  Figure  2A  and  2B  respectively.  In  Figure  2A,  AVT  increases  from  0.9V 

at  Vnr  =  0V  to  3V  at  VD„  —  12V.  The  increase  in  VT  is  due  to  the  widening 
ob  du  T,max  J 

of  the  depletion  width  by  the  backgate  bias.  In  Figure  2B,  the  decrease  in 
time  constant  is  due  to  optical  hole  emission  caused  by  the  external  IR 
radiation. 

EXPERIMENTAL  RESULTS 

Boron  doped  1-3  fl-cm  (100)  oriented  and  100-300  Cl-cm  (111)  oriented 
silicon  substrates  were  employed  in  this  study.  Gallium  doping  is  achieved 
by  using  Emulsitone*  spin-on  Galliumsilicafilm.  Pitting  of  the  substrate 
was  observed  after  Ga  diffusion.  N+  drain  and  source  regions  were  fabricated 
in  a  circular  geometry  using  Emulsitone  Phosphorsil icafilm.  Finished  devices 
were  then  mounted  on  the  cryo-tip  of  a  cryo-refrigeration  unit  (Air  Products 
Model  AC-2).  The  lowest  temperature  achievable  for  our  unit  is  20.4°K. 

A.  Thermal  Emission  Rate  of  the  Gallium  Center  in  the  Poole-Frenkel  Effect 
During  thermal  measurements  the  device  was  completely  shielded  to  prevent 
any  leakage  of  blackbody  radiation  and  stray  light.  Figure  3  shows  the  MOSFET 
characteristics  for  Ga  neutral  and  Ga  negative.  The  change  in  Ip<.  is  due  to 

thermal  hole  emission  from  the  Ga  centers. 

★ 
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Fig.  2  A.  Illustration  of  the  change  in  threshold  voltage  of  the  MOSFET  for  Ga  neutral 
and  Ga  negative  as  a  function  of  backage  gate  bias. 


The  hole  emission  rate  as  a  function  of  temperature  for  various  backgate 
bias  is  presented  in  Figure  4.  The  ionization  energy  of  the  Ga  center  can  be 
determined  from  the  relationship  [8], 

e£  =  AT3/2exp(-AE/kT)  (14) 

where  AE  is  the  observed  ionization  energy  of  the  Ga  centers,  k  is  the 
Boltzmann  constant,  T  is  the  temperature,  and  A  can  be  regarded  as  a  constant 
within  the  range  of  temperature  employed  to  determine  AE.  AE,  for  different 
backgate  bias,  is  given  in  the  inset  of  Figure  4.  At  a  backgate  bias  of  -12V, 

&E  =  63.3  meV,  which  is  close  to  the  published  value  of  65  meV  and  71  meV 
determined  by  the  Hall  and  optical  absorption  measurements  [13,14]. 

It  is  well  known  that  barrier  lowering  of  a  columbic  potential  due  to  an 
applied  e  field  gives  rise  to  field  enhanced  emission.  This  is  known  as  the 
Poole-Frenkel  effect  [15,16].  The  one  dimensional  Poole-Frenkel  effect  predicts 
an  enhancement  of  the  emission  rate  due  to  the  barrier  lowering. by  the  electric 
field  of  the  following  form  [15], 

ep(£)  =  ep(£  =  0)exp[B’/e/kT] 

=  AT3/2exp(-[AE(0)  -  B’v'el/kT)  (15) 

where  B’  =  (e3/TrK  eQ)1/2  =  2.2  x  10"4eV/cm7V,  AE(0)  is  the  Ga  ionization  energy 
at  zero  field,  Ks  the  dielectric  constant  of  silicon,  and  cQ  the  free  space 
permittivity. 

The  e  field  in  the  depletion  region  is  a  function  of  both  space  and  time. 
Since  e  is  not  a  constant  within  the  depletion  region,  an  average  e  field 
defined  as  e(t)  ■  4>S/W(t)  is  used.  The  surface  potential  <1>S  is  taken  at  L/2. 

<J>5  *  2<J>f  +  VQS  sat/2  -  VBG,  where  <J>f  is  the  Fermi  potential  in  the  bulk,  VDS>sat 


saturation  voltage  and  W  is  the  depletion  width.  For  case  (2),  because  Ng 
varies  from  at  t  =  0  to  Nj  at  t  =  »,  the  average  e  field  is  a  strong  function 
of  time.  In  view  of  Eq.  (15),  the  time  constant  observed  in  Eq.  (10)  is  weighted 
at  each  instant  of  time  by  the  e  field  at  that  instant.  The  data  were  measured 
between  t  =  O.lx  to  x.  To  show  the  effect  of  the  e  field,  we  have  used  the 

carrier  concentration  at  t  =  x/2  to  compute  the  e  field.  The  dependence  of 

mobility  on  the  e  field  also  introduces  an  additional  error  in  the  determination 
of  x. 

The  field  enhanced  hole  emission  rate  as  a  function  of  /E  for  various 
temperatures  is  presented  in  Figure  5.  S’  obtained  from  Eq.  (15)  is  given  in 
the  inset.  At  27°K,  the  measured  S’  is  approximately  a  factor  of  2  higher 
than  that  at  21 °K.  Since  the  e  field  in  the  depletion  region  is  both  a  function 
of  space  and  time,  the  measure  S’  is  at  best  an  average  value  which  does  not 
exclude  the  possibility  that  part  of  the  Ga  centers  beneath  the  inversion 
layer  could  have  tunneled  or  delocalized  under  high  field.  To  roughly  esti¬ 
mate  the  barrier  lowering  at  the  surface  using  the  theoretical  S’  =  2.2  x  ICf^eV 
(cm/V)^  at  Vgg  =  0  V  and  -12  V,  the  barrier  lowering  S ,y/£pea|<  1S  calculated 
to  be  79  meV  and  108  meV.  Compared  with  the  reported  gallium  ionization  energy 

of  63  meV  and  71  meV  using  the  Hall  and  optical  absorption  measurements,  the 

barrier  is  sufficiently  lowered  to  cause  tunnelling  and  delocalization  to 
occur.  The  actual  c  field  varies  from  a  maximum  below  the  inversion  to  zero 
at  the  depletion  edge.  This  fact  together  with  the  nonuniform  gallium  distri¬ 
bution  due  to  diffusion,  suggests  that  as  the  e  field  increases  only  those 
Ga  centers  at  the  smaller  value  of  the  e  field  (i.e.,  near  the  end  of  the 
depletion  region)  are  measured.  This  would  explain  the  increase  in  measured 
AE  as  the  e  field  increases.  Because  of  the  time  and  space  variation  of  the 
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electric  field,  our  measurements  are  not  accurate  enough  to  warrant  the  use 
of  the  three  dimensional  Poole-Frenkel  Effect  [17]. 


B.  Photoionization  of  the  Gallium  Center  in  Silicon 

When  the  MOSFET  device  structure  is  employed  as  an  infrared  detector, 
infrared  radiation  incident  on  the  gallium  centers  in  the  depletion  region 
causes  them  to  change  charge  state.  The  limiting  cases  where  all  of  the 
gallium  is  either  neutral  or  negative  has  been  shown  in  Figure  3.  The  ability 
of  an  impurity  center  to  capture  an  incident  photon  is  described  by  the  photo¬ 
ionization  cross  section.  In  order  to  more  fully  characterize  the  respon- 
sivity  of  the  IRFET  to  infrared  radiation,  one  should  determine  the  photo¬ 
ionization  cross  section,  o°(hu),  of  the  gallium  center,  which  is  related  to 
the  optical  emission  rate  through  [8]: 


ep  =  °°(h(i))$ 


(16) 


where  0  is  the  incident  photon  flux.  Since  the  observed  hole  emission  rate 
epbS  tlie  sum  t*ie  resP°nse  to  not  only  the  incident  signal,  but  also  the 
incident  background  flux  due  to  blackbody  radiation,  the  [18] 

b 


obs  _  o  ,  D. 
e„  =  e_  +  e 
P  P  P 


(17) 


where  e^  is  the  hole  emission  rate  due  to  the  background.  To  accurately 
determine  the  photoionization  cross  section,  the  intensity  of  the  incident 
monochromatic  signal  must  be  much  larger  than  the  total  background  flux  at 
all  wavelengths  to  which  the  gallium  center  is  responsive.  The  spectral  de¬ 
pendence  of  the  photoionization  cross  section  can  then  be  determined  by  varying 
the  incident  photon  energy  and  by  Eqs.  (1),  (2),  (3),  (16),  and  (17). 


A  theoretical  value  for  the  photoionization  cross  section  can  be  calculated 
using  the  Lucovsky  Model  [19,20] 

a°(hoj)  =  1  (eeff/eQ )2(167re2fi/3m*c) (AE)1/2(fia)  -  AE ) 3/2/ (fioa) 3  (18) 

where 

n  =  optical  index  of  refraction 

eeff/e0  =  effective  field  ratio 

m*  =  effective  mass  of  the  trapped  carrier 

AE  =  observed  ionization  energy 

c  =  speed  of  light. 

Or  in  terms  of  the  one  dimensional  Poole-Frenkel  effect: 

o°( hu)  -  l  (eeff/eo)2(16Tre2fi/3m*c)[AE(0)  -  B’v€]1/2[ti(JJ-(AE-e’^)]3/2/(1ia))3.  (19) 

As  pointed  out  in  the  last  section,  the  large  electric  field  in  the 
depletion  region  causes  a  reduction  of  the  observed  ionization  energy.  This 
will  undoubtedly  reduce  the  threshold  energy  of  the  photoionization  cross 
section.  However,  it  is  not  clear  whether  the  barrier  lowering  would  simply 
translate  the  entire  cross  section  to  the  lower  photon  energy  region. 

For  values  of  e  which  give  barrier  lowering  comparable  to  the  activation 
energy  AE(0),  tunnelling  or  delocalization  of  the  gallium  center  will  occur. 

This  means  that  the  surface  gallium  center  under  high  field  can  change  charge 
state  without  absorption  of  a  photon.  They  become  immediately  ionized  upon 
the  application  of  an  inversion  voltage  to  the  gate. 

A  300°K  background  radiates  large  photon  flux  beyond  4  microns  and  will 
discharge  the  Gallium  IRFET  in  tenths  of  milliseconds.  Our  experimental 
arrangement  for  measuring  the  photoionization  cross  section  consists  of  a 
globar  source  operating  at  1350°C  and  a  bandpass  filter  adjacent  to  the  device. 
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The  bandpass  filter  is  cooled  to  21 °K  and  does  not  constitute  a  source  of 
background  radiation  but  serves  to  block  off  the  background  radiation  outside 
the  bandpass  region.  This  increases  the  background  time  constant  to  the 
order  of  ten  milliseconds. 

Due  to  equipment  limitations  we  were  not  able  to  determine  the  absolute 
photoionization  cross  section  and  the  threshold  photon  energy.  Figure  6 
shows  the  relative  photoionization  cross  section  in  the  8  to  14  micron  wave¬ 
length  range  at  2.0. 4°K  for  three  electric  field  values.  The  dependence  of 
the  optical  emission  rate  on  the  average  electric  field  in  the  surface  de¬ 
pletion  region  is  shown  in  Figure  7.  By  comparing  Figure  5  &bd  Figure  7  we 
see  that  the  optical  emission  rate  exhibits  a  field  enhancement  similar  to  that 
observed  for  the  thermal  emission  rate. 

For  comparison,  measurement  of  Burnstein,  et.al.  [14],  theoretical  curves 
using  Eq.  (18)  are  also  plotted  in  Figure  6.  At  high  photon  energy  the 
measured  photoionization  cross  section  decreases  more  rapidly  than  the  theory 
predicts.  However,  the  Poole-Frenkel  barrier  lowering  does  seem  to  shift  the 
response  to  the  lower  photon  energy  region. 

OPERATION  AS  AN  INFRARED  DETECTOR 

In  previous  sections  an  investigation  has  been  made  of  the  character¬ 
istics  of  the  gallium  center  in  silicon.  When  operating  the  MOSFET  device 
structure  as  an  infrared  detector,  the  temperature  must  be  low  enough  so  that 
optical  emission  or  photoionization  is  the  dominant  process  as  opposed  to 
thermal  emission.  The  spectral  response  will  then  be  determined  by  spectral 
dependence  of  the  photoionization  cross  section  of  the  gallium  center. 

The  important  characteristics  describing  operation  of  the  device  as  an 
infrared  detector  are  then  the  quantum  efficiency,  responsivity  and  signal  to 


hi{j»6.  Relative  photoionization  cross  section  of  the  gallium  center  under  various 
electric  field  conditions.  Solid  curves  are  from  measurements.  Dashed  curves  are 
obtained  from  Lucovsky  Model.  Dash  and  dot  curve  was  measured  by  Bumstein  et.  al 
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When  the  device  is  operating,  the  quantum  efficiency  n  is  a  function 
of  time.  Ignoring  reflection  loss,  for  e°  »  e^,  n  is  given  by 

r  r 

n(t)  =  1  -  expC-o^je'^WCt)]  (20) 

W(t)  =  {2e0Kd>s/q[NA  +  Nj(l  -  e_t/T)]}1/2  (21) 

where  W  is  the  depletion  width  and  <j>s  is  the  surface  potential. 

At  time  t  =  0,  all  the  Ga  centers  are  neutral  and  able  to  emit  holes. 

At  this  time,  the  depletion  layer  width  determined  solely  by  the  number  of 
ionized  boron  atoms  is  the  widest.  According  to  Eqs.  (20)  and  (21)  the 
highest  quantum  efficiency  thus  occurs  at  t  =  0.  Upon  absorption  of  a  photon, 
the  neutral  Ga  emits  a  hole  and  becomes  negative.  As  the  photoionization 
process  progresses,  the  number  of  neutral  Ga  available  reduces  and  the 
depletion  width  decreases.  This  reduces  the  quantum  efficiency  rapidly. 
Specifically,  at  time  t,  such  that  Nj(l  -  e”t//T)  »  N^,  it  reduces  to 

n(t)  =  1  -  exp{-o°e"t/T[2c0Ks/q(l  -  e't/x)]1/2^}  .  (22) 

From  Eq.  (20)  or  Eq.  (22),  n(°°)  =  0,  i.e.,  all  the  Ga  centers  have  been 
ionized  and  no  further  photoionization  is  possible. 

Since  the  information  of  how  many  Ga  centers  have  been  photoionized  is 
stored  in  the  photocurrent  which  begins  at  t  =  0,  the  quantum  efficiency 
n(0)  can  be  taken  as  the  quantum  efficiency  of  the  IRFET.  From  Eq.  (20)  to 
achieve  high  n(0),  one  needs  to  have  high  gallium  concentration  and  wide 
depletion  width,  i.e.,  small  Ng  and  large  VgQ.  This  is  in  agreement  with  the 
results  obtained  in  Section  II-B  by  optimizing  the  photocurrent.  For 
o°  =  5  x  10’^cm2,  n(0)  as  a  function  of  Nj  for  various  Ng  and  VgQ  is  pre¬ 
sented  in  Figure  8.  In  this  figure,  the  effects  of  drain  to  source  voltage, 

breakdown  voltage  and  impurity  level  broadening  due  to  high  doping  levels 
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have  not  been  taken  into  account.  Employing  Nj  =  10  cm  ,  for  =  10  cm  , 
Vgg  =  20V,  an  initial  quantum  efficiency  of  23%  can  be  obtained. 

Since  the  amount  of  change  in  drain  to  source  current  depends  on  the 
length  of  integration  time,  for  a  given  $,  the  responsivity  can  be  defined  as 

r  s  change  in  photocurrent 
~  incident  photon  energy 


fiw$At 


I 


_  1  PS, max  "  Wt=V 


f)oj$At 


I 


(23) 

A/ Joule 


where  tj  is  the  integration  time,  ficu  is  the  photon  energy  and  A  is  the  gate 
area. 


To  provide  comparison  to  other  types  of  detectors,  a  responsivity  in 
terms  of  A/W  can  be  defined  as 


<24> 

where  we  have  assumed  that  the  photons  have  been  incident  on  the  device  until 

all  the  Ga  centers  have  been  ionized.  Of  course,  the  responsivity  would 

be  somewhat  unrealistic  since  for  a  given  integration  time  and  incident 

photon  flux,  the  Ga  centers  may  not  be  all  ionized.  Nonetheless,  in  the 

event  that  all  the  Ga  centers  will  be  ionized  or  that  tj  is  sufficiently  long 

such  that  Ine(O)  ■  Ine(tT)  -  AInc  „  ,  R0  does  indicate  the  IRFET's  ultimate 

05  os  I  OS, max  i 

performance.  Since  the  IRFET  can  be  designed  so  that  cutoff  occurs  before 
t  =  »,  which  could  be  considered  as  saturation  of  the  IRFET,  tj  could  be 
judiciously  chosen  for  a  particular  application. 
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It  has  been  previously  shown  that  the  maximum  attainable  current  signal 
to  noise  ratio  operating  under  shot  noise  limited  condition  is  [7] 

S/N  =  (N$A/2)1/2  (25) 

where  N<.  is  the  surface  charge  density  stored  by  the  Ga  center.  Shot  noise 
or  signal  noise  limited  operation  can  be  achieved  for  integration  period 
less  than  ^10  sec  or  for  short  integration  periods  such  that  the  1/f  noise 
is  smaller  than  the  shot  noise. 

High  responsivity  can  be  achieved  by  optimizing  the  photocurrent, 

Alps  max-  As  discussed  in  the  previous  section,  this  is  obtained  by  in¬ 
creasing  Vgg  and  Nj  until  e  ^  approaches  the  breakdown  field  value.  On 
the  other  hand,  too  large  an  e  field  will  cause  tunnelling  and  delocalization 
of  the  Ga  center  which  in  turn  reduces  Alp<.  max-  Heat  dissipation  and  device 
destruction  for  large  IpS,  should  also  be  considered. 

From  Eq.  (10),  after  AInc  has  been  optimized  it  can  be  kept  constant 

while  reducing  "A"  by  maintaining  the  same  W/L.  The  smallest  L  should  be 

large  enough  so  that  space-charge-limited  current  flow  will  not  occur.  From 

Eqs.  (22)  and  (23),  reducing  A  will  increase  R  but  reduce  S/N.  The  parameter 

A  may  be  reduced  to  the  extent  that  an  acceptable  S/N  ratio  is  obtained  as 

required  by  the  system  application.  The  parameters,  n,  Rj,  Rg,  and  S/N  for 

the  test  devices  are  given  in  Table  (1).  This  device  may  be  improved  by 

reducing  the  gate  area,  while  maintaining  W/L.  Parameters  for  an  improved 

design  are  also  listed  for  comparison.  In  the  table  a  flux  4>  =  3  x  16 
2 

photons/cm  *sec  at  X  *  14. 3y  is  assumed  incident  on  the  device.  For  the 
improved  design,  AVT  „„  increases  to  18.52  V.  A  current  of  69  mA  will 
undoubtedly  cause  device  destruction.  Although  the  responsivities  have 
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U  =  600  cm  vsec 
a  =  5  x  10~^8cm2 


XQ  =  14.3  ym 

4>  =  3  x  10^8  #/cm2sec 


PARAMETER 

PRESENT  DEVICE 

IMPROVED  DEVICE 

Xq(KA) 

1.75 

1.0 

W(ym) 

1.54  x  1 03 

200 

L(ym) 

2.285  x  102 

10 

NA(cm~3) 

5  x  10  13 

5  x  1013 

Nj(cm~3) 

3.4  x  1016 

6  x  1016 

iVT,max<V> 

2.66 

18.52 

AIDS,max^^ 

1.42 

68.85 

r/  (A/J) 

3.52  x  103 

1.68  x  108 

r2  (A/W) 

156.8 

8.4  x  106 

n(o)  (%) 

2.83 

6.84 

S/N 

1.3  x  105 

9.2  x  103 

W(0)  (ym) 

11.5 

23.6 

_ 

t  tj  is  taken  to  be  50  ms. 

Table  I.  Characteristics  of  the  current  devices  and  projected 
characteristics  of  a  better  small  area  device. 


.%  J 

■* .  •  %  «  "i 


V  • 


V  V 


increased  by  three  orders  of  magnitude  to  10  A/yw,  the  signal  to  noise  ratio 
has  dropped  an  order  of  magnitude  due  to  the  reduced  area.  The  increase  in 
quantum  efficiency  is  due  to  the  increase  in  backgate  bias. 

To  design  a  IRFET,  one  could  start  from  the  required  signal  to  noise 
ratio  dictated  by  the  particular  application.  From  the  expected  photo  flux  4>, 
the  desired  responsivity  and  integration  time,  the  parameters  Nj,  A,  W/L,  and 
Vg<.  are  then  determined.  Finally,  N^  can  be  determined  by  reset  time  require¬ 
ment. 

CONCLUSIONS 

In  conclusion,  operation  of  the  gallium-doped  MOSFET  as  an  infrared 
detector  in  8  to  14  ym  wavelength  region,  has  been  demonstrated.  Both  the 
characteristics  of  the  device  and  the  gallium  impurity  center  have  been 
described.  Measurements  of  the  thermal  emission  rate  of  the  gallium  center 
indicate  that  for  short  frame  times,  operation  at  20°K  should  be  acceptable, 
as  long  as  care  is  taken  not  to  subject  the  gallium  centers  to  large  electric 
fields  which  result  in  a  large  enhancement  of  thermal  emission  due  to  the 
Poole-Frenkel  effect.  The  time  constant  associated  with  thermal  emission 

4 

under  an  electric  field  of  2  x  10  V/cm  is  of  the  order  of  seconds  at  20°K. 

The  gallium  center  has  a  large  photoionization  cross  section  in  the 
far  infrared  region  and  if  exposed  to  radiation  from  a  300°K  scene  in  a 
large  aperture  or  small  f/number  system  will  completely  discharge  in  a  few 
milliseconds.  This  indicates  that  for  imaging  applications,  integration  and 
frame  times  in  the  tens  of  milliseconds  or  TV  rates  could  be  achieved.  The 
IRFET  has  been  found  to  possess  unique  characteristics  and  can  be  considered 
an  integrated  detector,  memory  element,  and  amplifier,  all  in  one  device;  as 
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opposed  to  CCD  scanned  photoconductor  arrays  where  each  component  is  a  separate 
element  of  device  [21-24].  As  such  then,  the  gallium  doped  infrared  sensing 
MOSFETs  might  be  particularly  suited  for  large  scale  integrated  infrared 
imaging  focal  planes  in  the  8  to  14  pm  wave  length  range. 
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NOISE  CHARACTERISTICS 

LIST  OF  SYMBOLS 

A 

2 

Area  of  individual  detector,  cm. 

AIP 

2 

Area  of  image  plane,  cm. 

BW,Av 

Bandwidth,  sec.”^  j 

c0 

_2 

Oxide  or  insulator  capacitance,  farads  cm.  ■( 

D* 

-1/2  -1  • 
Detectivity,  cm. sec  '  watt 

eo 

Electric  permittivity  of  free  space,  farad  cm.~^ 

•l 

Optical  emission  rate,  sec.~^  \ 

flOJ 

Signal  photon  energy,  joules  1 

(T2)1/2 
v  V 

Root  mean  square  (RMS)  noise  current,  amps. 

(T2  )V2 
UDSN; 

RMS  drain  to  source  noise  current,  amps.  { 

T 

Mean  or  average  current,  amps.  1 

JDS 

Drain  to  source  current,  amps.  if 

Ks 

Relative  dielectric  constant  of  silicon 

Ko 

Relative  dielectric  constant  of  oxide 

L 

Channel  length,  cm. 

n 

Number  of  point  in  image 

N 

-1  -2 

Background  photon  flux,  sec.  cm. 

na 

_3 

Acceptor  impurity  concentration,  cm. 

NEP 

Noise  equivalent  power,  watts 

Ns 

_o 

Surface  charge  density,  cm. 

ntt 

-3 

Impurity  concentration,  cm.  . 

q 

Electronic  charge,  coulombs 

(T*2)1'2 

RMS  variation  in  charge,  coulombs 

^CH 

Charge  in  channel,  coulombs.  I 
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Tj  Integration  time,  sec. 

Tp  Frame  time,  sec. 

Ts  Sample  time,  sec. 

S/N  Signal  to  noise  current  ratio 

VgG  Back  gate  or  substrate  bias,  volts 

VQS  Drain  to  source  voltage,  volts 

VGS  Gate  to  source  voltage,  volts 

VT  Threshold  voltage,  volts 

W  Gate  or  channel  width,  cm. 

Depletion  region  width,  cm. 
xQ  Oxide  thickness,  cm. 

A  Indicates  change  in  quantity  following  symbol 

n  Quantum  efficiency 

2  -1  -1 

\i  Mobility,  cm.  volt  sec. 

2 

a  Photoionization  cross  section,  cm. 

x  Time  constant,  sec. 

Fermi  potential,  volts. 

-1  -2 

4>  Signal  photon  flux,  sec.  cm 

Noise  equivalent  signal 
-1  -2 

photon  flux,  sec  cm 
INTRODUCTION 

The  ability  of  infrared  imaging  arrays  to  detect  signals  from  the  desired 
scene  is  limited  in  terrestial  applications  by  the  noise  resulting  from  the 
background  illumination.  Historically,  of  course,  the  simplest  technique  for 
infrared  imaging  has  been  to  mechanically  scan  the  scene  across  a  single 
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detector  area.  A,  where  the  bandwidth,  BW,  of  the  signal  amplification  system 
is  determined  by  the  number  of  points  in  the  scene  and  scan  rate.  This  leads 
to  the  definition  of  a  relative  figure  of  merit  for  comparing  detectors  as 
the  detectivity  or  0*  =  (A)^(BW)^ ^/NEP,  where  NEP  is  the  noise  equivalent 
power  [1].  In  the  case  of  the  newer  types  of  devices  as  for  instance,  charge- 
injection  diodes,  CID's,  [2,3],  charge  coupled  devices,  CCD's,  [3,4,5]  or 
infrared  sensing  field  effect  transistors,  IRFETs  [6,7,8,9,10]  there  is  no 
well  defined  single  bandwidth,  but  on  the  other  hand  two  bandwidths,  one 
associated  with  the  integration  period,  Tj,  and  one  associated  with  the  sample 
period,  T<.. 

The  first  objective  of  this  report  is  to  derive  the  noise  equivalent 
signal  power,  NEP,  for  an  array  of  IRFETs  employed  in  an  infrared  imaging 
array.  Since  we  are  dealing  first  with  the  theoretical  limit  on  the  operation 
of  infrared  imaging  detectors  due  to  background  radiation  it  will  be  assumed 
that  the  dominant  noise  mechnism  is  shot  noise  due  to  this  background.  Photon 
noise  will  not  be  considered  since  it  only  introduces  a  small  additional 
factor  similar  to  the  shot  noise  [1]  and  other  types  of  noise  as  1/f,  Johnson, 
and  generation  recombination  will  also  not  be  considered  at  first.  The  para¬ 
meters  for  the  imaging  application  are  the  area  of  the  scene  in  the  image 
plane,  AIp,  and  the  frame  time,  Tp,  background  photon  flux,  N,  number  of 
points  or  elements  in  the  scene,  n,  quantum  efficiency,  n»  and  signal  photon 
energy  fiw  .  Using  these  parameters  then,  the  theory  for  shot  noise  or  back- 
ground  limited  operation  of  the  IRFET  will  be  developed,  yielding  the  NEP  and 
signal  to  noise  ratio,  S/N. 

The  theoretical  expressions  for  shot  noise  limited  operation  of  the 
IRFET  will  then  be  compared  to  the  actual  operation  of  the  specific  example 


of  a  gold-doped  device  operated  in  the  near  infrared  range  at  78.8°K.  It 

will  be  shown  that  shot  noise  limited  operation  can  be  achieved  over  a  wide  ! 

range  of  operating  conditions.  In  reality,  however,  other  types  of  noise 

sources  are  also  present  in  the  MOSFET,  the  most  troublesome  being  1/f  noise. 

The  regions  of  shot  noise  and  1/f  noise  limited  operation  will  then  be  ' 

illustrated. 

SHOT  NOISE  THEORY  FOR  THE  IRFET 

The  infrared  sensing  MOSFET,  IRFET  [6,7,8,9,10]  will  integrate  the 
incident  flux  and  at  the  end  of  this  integration  period  there  will  be  a  change 
in  drain  to  source  current,  AIg<.(Tj),  in  the  linear  region  given  by  [6,8];  f 

I 

-AI0S(TI)  =  yC0(W/L)VDSAVTTIcN  (1) 

where  y  is  the  effective  surface  mobility,  Cq  the  oxide  capacitance,  W/L  the 
width  to  length  ratio  of  the  channel,  AVy  the  maximum  possible  change  in 
threshold  voltage,  a  the  photoionization  cross  section,  Tj  the  integration 
period,  and  N  the  background  flux.  However,  the  change  in  threshold  voltage 
is  due  to  a  change  of  charge  in  the  channel  device  and  there  will  be  associated 
with  this  some  uncertainty  in  this  change  in  charge  due  to  shot  noise.  As 
shown  in  Figure  1,  in  the  infrared  sensing  MOSFET,  the  emission  of  holes  from 

I 

impurity  centers  in  the  surface  space-charge  region  will  induce  a  change  of 
charge  in  the  channel  or  inversion  layer.  This  change  in  charge,  AQ^  can  be 
written 

15ch  ■  Tti 


where  I  is  the  average  source  current.  By  taking  into  account  the  displacement 


Fig.  1.  The  Infrared  Sensing  NDSFET  Showing  Source  and  Substrate  Current,  I,  due 
to  Emission  of  Holes  from  Impurity  Centers. 


current  or  change  in  width  of  the  surface  depletion  region  [11],  it  can  be 
shown; 


I  =  qtfNJTWdNWL/2 

(3) 

and  since  the  quantum  efficiency,  ,  is  just 

n  =  oNnwd 

(4) 

then 

T  =  qqWLN/2. 

(5) 

Since  this  is  a  time  dependent  process,  then  by  Carson's  Theorem  there  will 
be  a  root  mean  square  noise  current 

(^)1/2  =  (2qfAv)1/2  =  (q^WLN/Tj)1/2  (6) 

and  a  possible  error  in  the  channel  charge  at  the  end  of  this  integration 
period,  Tj  of 

(^)1/2  =  (q2nWLN/TI)1/2TI  .  (7) 

In  the  linear  region  of  operation 

!0S  =  yCo^W/L^VGS  "  VVDS  =  y(VDS/L2  ^ 

Qch  =  coulombs  (9) 

—2  1/2 

A  possible  error  of  charge  in  the  channel  (Q^)  at  the  end  of  the  integration 
period  will  result,  by  analogy  with  Eq.  8,  in  a  root  mean  square  noise  current 
in  the  sample  period  of 

<4n>1/2  ->MQjj>1/2v0S/L2  (10) 

(1D$N)1/2  ‘  M(q2nWLNTj)1/2V0s/L2 
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(11) 


due  to  the  signal  photon  flux  $  will 


The  signal  current,  Al0S(TI)  or  isigna-|» 
be  from  Eq.  (1) 

'signal  =  uCo<“/L>  VVtTl  <12> 

and  equating  the  signal  and  noise  to  determine  the  noise  equivalent  power 
yields  the  signal  photon  flux  which  gives  an  output  equal  to  the  noise 

$N  =  [l/(aCoAVT)][q2nN/(WLTI)]1/2  .  (13) 

However,  AVy,  can  be  expressed  in  terms  of  the  oxide  thickness  and  doping 
densities  [6,8] 

AVT  =  (Vs^o^2!^!  +  lVBGl)qNA/(,cseo)]1/2(NTT/2NA)  (14) 

then 

AVj  =  qWdNTT/(2Co)  (15) 

and  becomes 

4n  =  [4N/(WLnTj)]1/2  (16) 

and  the  noise  equivalent  power,  NEP  =  a>WL/or  a  single  detector 

NEP  =  [4NA/  (nTj )  ] 1  /2 .  (17) 

This  is  just  a  factor  of  (2)^2  larger  than  the  classical  formula  for  the 
shot  noise  limit  for  background  limited  operation  of  a  photodiode  [1]. 

In  an  imaging  application.  A,  the  area  of  each  detector  can  be  just 
A  *  Ajp/n,  and  the  integration  period  Tj  can  be  just  the  frame  time  Tp  and 
then  Tj  =  Tp.  Also  in  an  imaging  application  <j>  =  N.  The  NEP  for  a  single 
detector  or  element  in  the  scene  is  then 


NEP  =  fws[4NAIp/(nnTF)] 


(18) 
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This  is  a  factor  of  (n/2)  smaller  than  that  for  a  mechanically  scanned 
single  photodiode  detector  [12]. 

The  signal  to  noise  ratio  can  be  found  by  taking  the  ~atio  of  the  signal 
current  in  the  sample  period,  Eq.  12,  to  the  root  mean  square  noise  current 
in  the  sample  period,  Eq.  11,  to  give 

S/N  =  (nA$TI)1/2/2  (19) 

The  IRFET  will  saturate  if  the  product  of  $T^  is  too  large.  Saturation  will 
be  defined  to  occur  when  Tj  =  2t  =  2/($a)  =  2/e°,  and  usings  =  6  NTT  then 
the  maximum  attainable  signal  to  noise  ratio  is 

S"W  *  <W'V1/2-  <2°> 

The  product  of  is  just  the  change  in  charge  density  in  the  surface 

depletion  region,  if  we  let 

Ns  •  WdNn  (21) 


S'"*x  *  <NSA/2>1/2 


The  maximum  attainable  signal  to  noise  ratio  depends  only  on  the  stored  surface 
charge  density  and  area  of  the  device  under  shot  noise  or  background  limited 
operation.  It  is  independent  of  the  quantum  efficiency  integration  period, 

Tj  and  all  other  considerations.  Furthermore,  it  can  be  shown  that  a  similar 
expression  applies  to  CIO's  and  CCD's  as  well  so  that  all  surface  type  semi¬ 
conductor  devices  have  essentially  the  same  maximum  attainable  signal  to  noise 
ratio  since  they  all  have  comparable  surface  charge  storage  densities  which 
are  limited  by  the  dielectric  strength  of  the  gate  insulators  or  avalanche 
breakdown  in  the  substrate. 


EXPERIMENTAL  APPARATUS 


The  device  used  in  the  measurements  was  made  on  (100)  orientation  silicon 
p-type  substrate  of  l-2ft‘cm  resistivity  and  doped  with  a  gold  concentration 
of  2  x  1015  cm-^.  The  gate  oxide  of  the  MOSFET  is  4000  A  thick  and  the 
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circular  channel  has  an  8.17  W/L  ratio  providing  an  active  area  of  2.5  x  10  cm 
[7,8].  By  mounting  the  device  on  a  T0-5  header  fitted  in  a  tubular  reflective 
cavity,  both  frontside  and  edge  illumination  was  possible.  The  details  of 
operation  of  this  specific  example  of  a  gold  doped  silicon  IRFET  employing 
the  response  of  the  donor  level  has  been  treated  in  detail  in  the  literature  [8] 
In  this  case,  x,  the  time  constant  associated  with  the  change  in  charge  state 
of  the  gold  donor  level  will  be  designated  as  x^,  which  is  the  reciprocal  of 
the  product  of  the  incident  photon  flux  and  photoionization  cross  section  [8]. 

The  experimental  arrangement  and  the  equipment  used  in  these  measurements 
are  shown  in  Fig.  2.  The  upper-corner  frequency  of  the  Hewlett-Packard  DC 
MicroVolt-Arrmeter  425A  is  l-Hz(-3db),  and  the  lower-corner  frequency  of  the 
Krohn-Hite  3750  Filter  (high  pass)  is  adjustable  from  0.02-Hz;  therefore, 
the  two  instruments  provide  for  a  variable  noise  bandwidth.  A  Bausch  and 
Lomb  tungsten  light  source  powered  from  a  variac  permitted  various  photon 
flux  intensities.  The  output  signal  current  and  the  accompanying  noise  was 
measured  with  the  IRFET  running  continuously  and  recorded  by  a  Houston  Instru¬ 
ment  Omnigraphic  2000  X-Y  Recorder.  A  1-uF  DC-blocking  capacitor  was  attached 
to  the  input  of  the  X-Y  recorder  (1-MH  input  impedance)  in  order  to  eliminate 
any  drift  in  the  high-pass  filter/volt  meter  output.  This  circuit  attenuated 
all  frequencies  below  0.159-Hz,  thus  limiting  the  magnitude  of  the  1/f  type 
noise.  The  incident  photon  flux  was  determined  by  noting  the  time  tp  required 
for  lAlp^l  to  increase  from  0.1  to  0.9  of  its  final  value  and  this  is  shown 
in  Figure  3. 


Experimental  Apparatus  used  in  the  Noise  Measurements 


EXPERIMENTAL  RESULTS 


The  overall  system  noise  with  the  high-pass  filter  input  shorted  was 

about  8mV  „  at  the  voltmeter  output.  The  IRFET  was  connected  in  a  common 
rms 

source  configuration  with  a  precision  1  kn  load  resistor. 

A.  Interface-State  Noise 

The  flow  of  electrons  from  one  site  to  another  in  the  surface  channel  of 

a  MOSFET  is  a  random  process  because  of  trapping  and  emission  of  carriers  by 

interface  states  [13,14].  Slow  interface  states  (oxide  states)  produce 

appreciable  noise  in  IRFETs  during  long  integration  periods.  Whenever  the 

device  is  active,  both  fast  and  slow  interface-state  noises  are  present  as 

a  dark  current  with  1/f-type  frequency  dependence  [13,14], 

This  dark-current  noise  was  observed  in  different  bandwidths  for  170  sec 

by  varying  the  lower  corner  frequency  of  the  high-pass  filter  from  0.16  to 

4  Hz  (10  db  gain  after  1  Hz)  and  recording  the  RMS  noise  voltage  across  the 

1  kfl  load  resistor  versus  time.  By  taking  the  difference  in  noise  current 

between  two  observations  differing  by  Af  =  0.2  Hz  and  dividing  by  0.2  Hz, 

the  dark-current  noise  was  specified  for  a  given  center  frequency  in  a  0.2  Hz 

bandwidth.  The  results  of  this  series  of  measurements  are  shown  in  Figure  4 

(note  that  1/f  =  Tj).  Fitting  the  experimental  points  to  a  straight  line 

(method  of  least  squares)  indicates  that  the  dark-current  noise  in  a  cooled 
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IRFET  varies  approximately  as  (1/f)  '  . 

B.  Shot  Noise 

By  illuminating  the  IRFET  with  infrared  light,  the  shot  (or  quantum) 
noise  associated  with  the  optically  generated  change  in  the  channel -charge 
could  be  measured.  A  typical  result  comparing  the  measured  dark-noise  and 
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shot-noise  currents  in  a  bandwidth  Av  =  0.2  Hz  is  shown  in  Figure  5.  During 
illumination  the  noise  current  is  seen  to  eventually  approach  the  dark  noise 
level  as  the  length  of  the  integration  time  increases  and  the  gold  centers 
discharge.  In  order  to  observe  the  initial  shot  noise  current  immediately 
after  the  IRFET  was  illuminated,  the  time  sacle  can  be  expanded. 

The  expected  shot  noise  current  can  be  calculated  by  using  Eqs.  11  and 
4,  and  the  expression  for  the  maximum  possible  change  in  drain  current  [8] 

AlJs  -  UC0(W/L)4VTVDS  (23) 

which  becomes  after  using  Eq.  15 

«DS  ■  »C0WL)tqWdNTTVDS/(2C0)]  (24) 

and  then  by  using  Eq.  11  and  8 

^DSN  ^/2  =  £2qAlDS  VGS  "  VA)11/2- 

In  reality,  this  expression  has  been  developed  assuming  that  the  IRFET  is 
operating  in  the  linear  region  and  a  linear  approximation  has  been  applied 
to  describe  the  time  response.  However,  it  will  be  used  to  calculate  and 
show  that  the  theoretical  shot  noise  is  of  the  same  order  as  that  observed 
in  operation  of  the  device  even  though  the  above  conditions  might  not  be 
rigidly  adhered  to.  For  Figure  5,  1^  is  1.8  mA,  AIp<-  is  0.4  mA,  Tj  a 
100  seconds,  rn  is  100  seconds,  C  can  be  calculated  from  the  oxide  thickness 
of  4000  A,  the  area  A  is  2.5  x  10"^  cm^,  and  (Vgs  -  Vy)  is  12.  5  volts.  This 
yields  a  root  mean  square  noise  current  of  30  nA  or  a  peak  to  peak  noise 
current  of  85  nA,  when  substitutions  are  made  into  Eq.  25.  The  actual  shot 
noise  observed  in  the  first  few  seconds  in  Figure  5  after  illumination  is 
less  than  the  order  300  nA  peak  to  peak.  Note  that  the  time  scale  over 


which  the  current  fluctuations  are  being  plotted  in  Figure  5  is  ten  sec.  per 
division.  This  rules  out  the  use  of  a  meter  to  determine  the  RMS  value  of 
the  noise  signal.  The  measured  shot  noise  is  expected  to  be  somewhat  larger 
than  the  calculated  value  due  to  the  fact  that  the  actual  device  is  being 
operated  in  the  saturation  region  rather  than  the  linear  region  for  which 
the  theory  has  been  developed  to  simplify  the  analysis. 

The  quantity  (T^  )^,  given  by  Eq.  25  represents  the  uncertainty 
which  occurs  in  the  drain  current  at  the  end  of  the  integration  period,  or 
at  the  beginning  of  a  sample  period  as  might  occur  in  an  application  to  read 
the  conductivity  state  [6].  The  fluctuations  in  Figure  5  during  illumination 
are  representative  of  this  uncertainty.  Note  that  the  input  circuit  of  the 
X-Y  recorder  limits  1/f  type  circuit  noise  to  40  nA  peak  to  peak  or  approxi¬ 
mately  12  nA  RMS,  due  to  the  cut-off  frequency  of  0.159  Hz  or  time  constant 
of  6.3  sec.  Thus,  shot  noise  could  still  be  observed  for  long  integration 
periods.  The  expected  actual  noise  (1/f- type  dominant)  as  would  be  seen  by 
a  static-change  (DC)  meter  is  much  larger  as  will  be  shown  in  Figure  6. 

Note  that 

td  =  1/No  (26) 

so  that  in  Eq.  25  then,  (T^-N  )^  will  increase  with  the  incident  photon 
flux.  If  the  photon  flux  is  kept  constant,  then  the  noise  current  increases 
with  integration  time  in  Eq.  25.  The  apparent  decrease  in  shot  noise  is  due 
to  the  actual  non-linear  time  response  of  the  device  and  discharge  of  an 
appreciable  number  of  gold  centers  for  times  of  the  order  ip.  A  signal -to- sKot 
noise  ratio  can  be  calculated  using  Eq.  25  and  the  fact  that 
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to  yield 


S/N  =  CCo(VGS  '  vj)/2cl)(AlDS/IDS)(TI/TD)A]1/2-  (28) 

The  signal  to  noise  ratio  Increases  with  both  the  ratio  (Tj/Tq)  and  the  area 

A  of  the  device.  (Note  the  similarity  to  Eq.  22  for  Tj  =  2tq.) 

The  expected  drain-to-source  RMS  noise  current  has  been  plotted  in 

Figure  6  for  different  integration  times  as  a  function  of  incident  photon 

flux.  The  shot  noise  for  a  100  second  integration  period  is  shown  as  a 
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dotted  line  and  the  value  of  300/2  v?  nA  at  a  flux  of  6.67  x  10  photons/cm 
sec.  taken  from  Figure  5.  The  expected  1/f  noise  in  this  case  is,  however, 
much  larger  as  obtained  from  extrapolating  Figure  4  as  1800  nA  and  is  shown 
as  a  solid  line  in  Figure  6.  For  shorter  integration  periods  the  1/f  noise 
will  be  much  smaller  and  the  expected  noise  characteristics  for  integration 
periods  of  1  and  10  seconds  will  be  shot  noise  limited  at  the  higher  photon 
fluxes. 

Figure  6  has  been  divided  into  three  distinctly  different  regions  which 
are  representative  of  the  different  operating  conditions.  In  the  saturated 
region  the  photon  flux  Is  too  large  and  the  integration  period  too  long  to 
allow  operation  of  the  device  as  a  photon  detector,  all  the  impurity  centers 
discharge  before  the  end  of  the  Integration  period.  In  the  1/f  noise  limited 
region,  which  corresponds  to  low  photon  fluxes  and  long  integration  period, 
1/f  rather  than  shot  noise  Is  the  limiting  factor.  In  the  shot  noise  or 
background  limited  region,  which  corresponds  to  relatively  higher  photon 
fluxes  and  shorter  integration  periods,  shot  noise  due  to  the  incident  photon 
flux  Is  the  dominant  noise  mechanism.  It  Is  shown  that  shot  noise  or  back¬ 
ground  limited  operation  of  the  IRFET  can  be  achieved  for  Integration  periods 
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of  less  than  about  10  seconds  or  with  photon  fluxes  above  about  10  cm 


T 


and  similar  characteristic  curves  would  apply  to  other  individual  devices 
and  no  special  effort  has  been  made  to  minimize  the  surface  state  density 
and  1/f  noise  in  these  devices. 

APPLICATION  CONSIDERATIONS 

The  photon  flux  incident  on  a  particular  detector  should  be  anticipated 
to  vary  by  orders  of  magnitude  depending  upon  the  wavelength  and  type  of 
application.  In  the  case  of  gold-doped  devices  [8],  it  might  be  anticipated 
that  in  the  near  infrared  photon  fluxes  may  well  be  low  unless  the  detector 
is  being  employed  to  localize  a  very  hot  source  and  the  1/f  noise  limit  may 
be  an  important  consideration.  For  gallium  doped  devices, on  the  other  hand 
[10],  employed  in  terrestial  applications,  the  flux  incident  on  the  detector 

is  N  =  N300^Ac^^4f2^  N30o(Xc )  the  blackbody  flux  from  a  300°K  source 

integrated  out  to  the  longest  wavelength  of  response  of  the  impurity  center, 
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17  pm  in  this  case  [1].  With  F  =  1.5  optics  N  =  2  x  10  /cm  sec  and  the 

time  constant  of  discharge  of  the  gallium  centers  would  be  x  =  l/(oN)  or 

approximately  16  milliseconds  if  an  average  photoionization  cross  section  of 
-16  2 

a  =  3  x  10  /cm  is  used.  In  the  case  of  gallium  doped  devices  in  terrestial 
applications,  integration  periods  or  frame  times  might  be  anticipated  to  be 
short  and  the  shot  noise  or  background  limit  apply.  In  this  case,  the  I/f 
noise  limit  is  of  little  concern.  Indium-doped  devices  [9,10,15]  are 
intermediate  to  the  above  cases.  In  terrestial  applications  N^gg^um)  = 

2  x  1017  and  using  an  average  photoionization  cross  section  of  1.5  x  10’^/cm^ 
and  F  =  1.5  yields  a  time  constant  of  300  milliseconds. 

Different  considerations  would  of  course  apply  in  extra-terrestial  appli¬ 
cations  and  low  flux  densities  may  be  incident  at  all  wavelengths,  in  which 
case  long  integration  periods  might  be  employed  and  the  1/f  noise  limit  applied. 


CONCLUSIONS 


It  has  been  shown  theoretically  that  under  shot  noise  or  background 
limited  operation  that  the  maximum  attainable  signal  to  noise  ratio  in  the 
infrared  sensing  MOSFET  (IRFET)  is  determined  only  by  the  area  of  the 
detector  and  the  surface  charge  storage  density.  Equations  describing  the 
noise  characteristics  of  this  new  type  of  integrating  infrared  detector  have 
been  developed  and  presented,  assuming  operation  of  the  device  in  the  linear 
region  and  employing  a  linear  approximation  to  the  time  response  of  the 
device. 

Noise  measurements  on  the  specific  example  of  a  gold  doped  device  employ¬ 
ing  response  of  the  gold  donor  level  have  identified  the  shot  noise  associated 
with  the  incident  infrared  signal  and  the  1/f  noise  inherent  in  the  device 
itself.  The  regions  of  shot  noise  or  background  limited  operation  and  1/f 
noise  limited  operation  have  been  delineated. 

Although  measurements  have  been  made  only  on  gold  doped  devices,  similar 
results  might  be  anticipated  for  both  the  indium  and  gallium  doped  infrared 
sensing  MOSFETs.  These  results  indicate  that  for  shorter  integration  periods 
of  around  10  seconds  or  less,  shot  noise  or  background  limited  operation  of 
the  IRFET  should  be  possible  in  infrared  detector  and  imaging  applications. 
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UNIFORMITY  OF  DOPINGS 


A  number  of  indium-  and  gallium-doped  wafers  were  oxidized  and 
detailed  measurements  made  on  MOS  capacitors  to  determine  the  uniformity 
of  dopings  being  acheived  in  the  wafers  since  the  device  reponse  depends 
directly  on  this  doping.  The  IRFET  is,  however,  not  like  photoconductors 
in  that  in  the  IRFET  the  responsivity  does  not  depend  upon  carrier  lifetimes 
and  concentrations  of  residual  impurities  in  the  wafers.  As  such  then  it 
is  only  necessary  to  determine  the  uniformity  of  the  main  or  majority  dop¬ 
ant  in  the  wafer,  either  indium  or  gallium. 

Fig.  1  shows  the  spatial  variation  of  doping  across  a  gallium-doped 
wafer.  Fig.  2  presents  histograms  of  the  variation  in  doping  concentrations 
on  both  indium-  and  gallium-doped  wafers.  The  results  show  that  a  variation 
in  doping  by  a  factor  of  "two"  is  observed  across  both  types  of  wafers. 

The  indium-doped  wafers  were  obtained  from  General  Diode,  Framingham, 
Mass,  as  an  off  the  self  stock  item.  The  gallium  doped  wafers  were  received 
from  a  DOD  contractor  which  had  some  left  over  from  a  prior  project.  In 
both  cases  no  special  care  was  taken  to  insure  a  uniform  radial  distribution 
of  dopings  on  the  wafer.  These  results  correspond  to  the  previous  integrated 
circuit's  experience  of  the  Principal  Investigator  were  is  was  found  in 
commercial  boron-doped  wafers  used  in  MOSFET  integrated  circuits  that  there 
was  in  general  a  variation  of  dopings  by  a  factor  of  "two"  in  wafers  of  the 
same  lot  and  radial  variations  of  a  factor  of  "two". 

In  addition  to  doping  variations  the  MOSFET  device  structure  is  also 
sensitive  to  variations  In  threshold  voltage.  Experience  has  shown  that  even 
on  the  same  chip  threshold  voltages  would  vary  by  0.1  Volt.  As  such  then  the 
level  of  uniformity  of  responsivity  which  can  probably  be  acheived  with  some 
care  for  MOSFET  detectors  on  the  same  chip  is  of  the  order  10.0%.  Clearly  in 
most  applications  signal  processing  will  be  required;  as  was  proposed  origin¬ 
ally  this  would  involve  the  use  of  auxiliary  memory  planes. 
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APPLICATION  CONSIDERATIONS 


APPLICATION  TECHNIQUES 

As  was  described  at  previous  points  in  the  text  the  IRFET  detector 
would  be  employed  in  an  X-Y  addressed  array  as  shown  in  Fig.  1.  This  imaging 
array  essentially  constitutes  a  read-only  memory  plane  where  the  conductance 
of  each  device  is  a  function  of  the  infrared  flux  incident  on  that  partic¬ 
ular  location.  Associated  with  this  imaging  memory  plane  would  be  off-chip 
logical  control  elements  and  auxiliary  memory  planes  for  signal  processing  as 
shown  in  Fig.  2. 

The  conductance  of  each  element  in  the  imaging  array  would  be  read  on 
an  analog  basis  both  before  and  after  illumination  and  converted  to  a  digital 
format  by  the  A/D  converter  shown  in  Fig.  2.  In  doing  so,  note  that  since  the 
IRFET  is  a  static  or  D.C.  memory  element  the  read  operation  can  be  suppressed 
until  the  address  location  becomes  stable.  Following  this  the  length  of  the 
read  operation,  or  sample  time,  is  determined  by  the  precision  which  is  desired 
or  required  in  reading  the  conductance  of  the  element.  In  theory  these  sample 
periods  could  be  made  very  long  allthough  at  the  expense  of  long  frame  times. 

The  random-access  memories  (RAM's)  in  Fig.  2  would  be  used  to  store 
the  conductance  of  each  imaging  element  or  IRFET  before  and  after  illumination. 
The  programable  read-only  memory  (p-ROM)  would  be  used  to  compensate  for  any 
non-uniformities  in  the  responsivity  of  elements  in  the  array.  The  net  result 
would  be  to  allow  a  very  high  degree  of  resolution  between  different  flux  inten¬ 
sities  at  different  imaging  elements  to  be  acheived. 

This  signal  processing  in  some  senses  is  analogous  to  the  "time-delay  and 
integration  (TDI)"  used  in  CCD  scanned  photoconductors.  The  difference  here  is 
due  to  the  fact  that  the  IRFET  works  on  a  static  as  opposed  to  dynamic  memory 
concept.  Note  also  that  the  IRFET  array  could  acheive  a  very  high  resolution 
in  a  full  staring  mode.  The  IRFET  array  and  signal  processing  technique  does 
not  require  mechanical  motion  of  either  the  target  or  the  array,  as  opposed 

to  TDI  techniques  which  require  motion. 
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COMPARISON  TO  CCD  SCANNED  PHOTOCONDUCTOR  ARRAYS 

Table  I.  presents  a  comparison  of  the  operation  characteristics 
of  a  32  bit  CCD  scanned  array  of  gallium-doped  photoconductors[7,8]  to 
the  projected  characteristics  of  an  IRFET  array  based  on  our  current 
devices[l-6]  and  based  on  better  smaller  area  devices.  Such  a  comparison 
is  difficult  to  make  since  the  two  different  types  of  detectors  have  such 
different  operating  principles  and  it  must  be  kept  in  mind  in  doing  so 
how  each  detector  works. 

While  this  Table  summarizes  a  number  of  points  from  previous  sections 
a  number  of  points  are  worthwhile  considering  again  in  detail.  Note  the 
most  striking  feature  is  that  the  IRFET  detectors  have  a  much  higher  respon- 
sivity,  and  which  was  pointed  out  previously  does  not  depend  upon  carrier 
lifetimes  and  residual  impurity  concentrations.  Since  the  IRFET  employs 
only  the  gallium  centers  within  about  one  mil  from  the  surface  the  quantum 
efficiency  is  on  the  other  hand  somewhat  lower. 

The  theoretically  maximum  attainable  signal  to  noise  ratio  in  both 
types  of  systems  is  essentially  just  the  square  root  of  the  total  number  of 
stored  electrons.  In  the  IRFET  detector  the  detector  and  integration  element 
are  one  and  the  same  and  have  the  same  area.  In  the  CCD  system  the  storage 
is  usually  smaller  than  the  detector.  Consequently  the  IRFET  can  acheive 
at  least  the  same  if  not  a  significantly  larger  signal  to  noise  ratio. 

The  bottom  lines  in  the  Table  emphasis  the  difference  in  the  read-out 
techniques  and  the  signal  at  the  output  of  the  system.  In  the  CCD  system  the 
output  is  an  integrated  charge  in  coulombs;  in  the  IRFET  system  it  would  be  a 
current  which  can  be  read  for  a  Indefinitely  long  time  period  on  a  non-destr¬ 
uctive  basis.  In  the  IRFET  the  integrated  charge  does  not  flow  in  the  external 
circuit  but  rather  remains  stored  in  the  depletion  region  of  the  MOSFET  until 
the  reset  operation.  The  Table  also  indicates  that  acceptable  frame  rates  could 
be  acheived  in  most  applications  using  the  IRFET  detector  array. 
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A  COMPARISON  OF 

GALLIUM-DOPED  SILICON  INFRARED  DETECTORS 
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Table  I.  A  comparison  of  the  characteristics  of  CCD  scanned  photoconductors 

to  infrared  sensing  MOSFET's. 
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CONCLUSIONS 


This  report  provides  a  detailed  description  of  the  characteristics 
of  gold,  indium,  and  gallium  doped  MOSFET  device  structures  as  infrared 
detectors.  In  doing  so  it  also  provides  a  detailed  description  of  the 
characteristics  of  these  three  impurities  in  silicon.  In  the  case  of  gold, 
the  results  obtained  could  be  compared  to  those  obtained  previously  by  this 
and  other  authors.  In  the  case  of  indium  and  gallium  the  work  was  much  more 
exploratory  and  this  work  constitutes  to  the  author's  best  knowledge  the 
first  direct  observation  of  the  thermal  and  optical  emission  characteristics 
of  these  centers. 

Work  was  begun  on  gold-doped  devices  first  because  of  the  author's 
prior  experience  with  this  impurity  center  and  due  to  the  fact  that  the 
measurements  could  be  easily  done  above  77  °K.  At  the  beginning  of  the  project 
facilities  existed  only  for  measurements  above  liquid  nitrogen  temperatures 
and  the  laboratory  facility  at  the  University  of  Arkansas  for  processing 
MOSFET  devices  was  just  being  set-up.  This  work,  on  gold-doped  devices,  served 
to  demonstrate  that  the  device  did  if  fact  work  according  to  the  original 
design  concepts  which  had  been  previously  published  in  August  1974.  The 
gold-doped  devices,  however,  had  a  low  quantum  efficiency;  primarily  since 
it  was  desired  not  to  overcompensate  the  substrate  and  make  it  highly  resist¬ 
ive  at  room  temperature. 

Following  this  work  was  initiated  on  indium  and  gallium  doped  devices. 

By  this  time  sufficient  experience  had  been  obtained  in  using  the  MOSFET  device 
structure  as  an  infrared  detector  to  allow  a  design  and  acquisition  of  the  low 
temperature  and  long  wavelength  infrared  systems  required  in  characterizing  the 
device.  Classical  chopped  radiation  systems  are  not  approplate  for  an  Integrating 
infrared  detector  and  many  of  our  arrangements  are  rather  unique.  In  addition 
*' —  Anna/t  4nH4nm  waferc  had  been  obtained  and  a  diffusion  source  for 


neither  indium  nor  gallium  doped  wafers  were  readily  available,  even  from 
D00  contractors  who  were  also  just  starting  to  become  actively  engaged  in 
work  on  indium  and  gallium  doped  silicon  detectors. 

The  subsequent  work  on  indium  and  gallium  doped  detectors  proceeded 
in  a  relatively  straightforward  manner.  A  very  significant  difference  was 
however  to  occurr  in  later  work  with  the  realization  that  lightly  boron 
doped  substrates  could  be  employed  to  make  a  very  large  improvement  in  the 
quantum  efficiency  of  the  detectors.  One  objection  to  use  of  the  IRFET 
detector  had  been  the  low  quantum  efficiency  reported  in  the  initial  work 
and  in  the  original  design  proposal.  Using  devices  with  a  large  indium  or 
gallium  concentration  and  a  low  boron  concentration  necessitated  reworking 
all  equations  describing  device  operation  and  the  development  of  different 
analytical  techniques.  The  work  presented  here  is  in  the  historical  order, 
where  the  introduction  serves  to  outline  all  the  different  possibilities  for 
doping  ratios  and  device  reponse  characteristics. 

If  one  Is  most  interested  in  impurity  center  characteristics,  for  the 
measurements  at  least,  It  Is  more  convenient  to  have  a  large  boron  concentr¬ 
ation  and  low  indium  and  gallium  concentration.  Such  crystals  would  however 
be  hard  to  grow. 

This  report  we  believe  demonstrates  that  the  IRFET  device  structure 
can  in  fact  be  used  as  a  detector  in  the  near,  middle,  and  far  infrared  wave¬ 
length  regions.  It's  very  unique  characteristics  might  make  it  particularly 
advantageous  in  some  applications. 
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APPENDIX  A:  Application  of  the  MOSFET  Device  Structure  in  Characterizing 
Impurity  Centers _ 

Identification  of  impurity  center  concentrations,  emission  rates  and 
capture  coefficients  using  pulsed  capacitance  techniques  on  p-n  junctions, 
Schottky  barrier  junctions,  and  MIS  capacitors  are  low  temperature  measure¬ 
ments.  As  shown  in  Fig.  1,  these  structures  are  two  terminal  devices  which 
are  in  series  with  the  substrate  resistance. 

In  the  case  of  indium  and  gallium  doped  devices,  the  required  temperature 
of  operation  are  in  the  range  of  50°K  to  20°K  and  a  very  significant  deioni¬ 
zation  will  occur  in  the  substrate.  This  results  in  a  very  high  series 
impedance,  Rsubstrate’  ma^es  Pulsed  high  frequency  capacitance  measure¬ 

ments  virtually  impossible.  While  such  capacitance  measurements  can  be  made 
on  a  double  indium-boron  doped  substrate,  in  the  case  of  gallium,  however, 
even  a  double  gallium-boron  doped  substrate  will  not  work  since  even  deioni¬ 
zation  of  the  boron  becomes  appreciable  at  the  required  temperatures  of  opera¬ 
tion  resulting  in  an  excessive  substrate  impedance. 

The  MOSFET  structure  on  the  other  hand  is  a  3  terminal  measurement  (ex¬ 
clusive  of  the  gate)  where  the  surface  conductivity  between  source  and  drain 
is  measured.  The  substrate  contact  serves  only  to  establish  the  substrate 
potential  and  at  low  temperatures  only  an  unmeasureably  small  leakage  current 
is  drawn  through  this  substrate  contact.  As  a  consequence,  measurements  can 
still  be  easily  made  on  the  MOSFET  structure  when  the  substrate  is  deionized 
and  high  Impedance  ('v  1  megohm  or  more)  since  almost  no  potential  drop  occurs 
across  the  substrate.  The  MOSFET  structure  is  thus  ideally  suited  for 
measurements  on  single  indium-doped  or  gallium-doped  materials  with  only 
residual  impurity  dopings  of  shallow  level  acceptors  or  donors  as  boron  or 
phosphorus,  such  as  ones  used  in  extrinsic  silicon  photoconductive  detectors. 

In  fact,  the  concentration  of  such  residual  impurities  can  be  easily  determined 


Substrate  Contact 
for  MOSFET 


Figure  1  Devices  Used  for  the  Characterization  of  Impurities  in  Silicon; 

In  the  upper  three,  p-n  junction,  Schottky  Barrier  Junction,  and 
MIS  capacitor  the  capacitance  of  the  device  is  in  series  with  the 
substrate  resistance  which  becomes  much  larger  than  the  impedance 
of  the  depletion  region  capacitance  at  low  temperatures  on  a 
deionized  Substrate,  raking  capacitance  measurements  very  diffici 

In  the  MOSFET  structure  essentially  no  current  is  drawn  through  ' 
substrate  resistance  and  substrate  contact. 


as  will  be  discussed. 


A.  Illustration  of  Capability  of  MOSFET 

The  results  reviewed  in  this  section  were  obtained  on  a  single  indium 
doped  silicon  substrate  such  as  that  which  might  well  be  employed  in  a 
extrinsic  silicon  photoconductive  detector  for  the  middle,  3  to  5  micrometer, 
infrared  wavelength  range.  Specifically,  results  are  shown  for  a  character¬ 
ization  of  the  indium  center  in  silicon  and  since  these  results  have  been 
treated  to  some  extent  in  the  literature,  a  detailed  description  will  not 
be  given.  Sub-section  (iv),  describes  the  proposed  technique  which  could  be 
employed  to  determine  the  carrier  capture  coefficient  or  cross  section  an 
impurity  center  using  the  MOSFET  structure. 

(i)  Impurity  Concentrations 

Impurity  concentrations  are  easily  determined  employing  the  MOSFET 
structure  by  measuring  the  dependence  of  threshold  voltage  on  back-gate 
basis.  Fig.  4  on  page  55  of  this  report  is  an  illustration  of  this  technique 
on  an  indium  doped  substrate.  This  particular  sample  has  an  indium  concen- 

1C  1 

tration  of  8.92  x  10  /cm  and  a  net  residual  shallow  level  acceptor  con¬ 
centration  of  4.95  x  lO^/cm'*. 

(ii)  Thermal  Emission  Rates 

Fig.  6  on  page  59  of  this  report  shows  a  determination  of  the  thermal 
emission  rate  by  measuring  the  temperature  dependence  of  the  time  constant  of 


the  dark  decay  between  the  two  limiting  conductivity  states.  The  slope  of  the 
line,  0.167eV,  is  appropriate  to  the  thermal  activation  energy  of  the  indium 
center  in  silicon  and  provides  a  unique  determination  and  identification  of 
the  indium  center. 

The  emission  rate  of  shallow  level  centers  determines  the  rate  of  ioni¬ 
zation  of  the  centers  as  a  function  of  temperature  in  photoconductive  detec¬ 
tors  and  in  general  have  not  been  previously  determined,  although  obviously 
very  important. 

This  sample  was  very  carefully  checked  for  other  levels,  as  for  instance 
the  O.lleV  level  reported  by  Hughes  in  indium  doped  samples,  however  no  such 
other  levels  were  found  in  our  samples.  The  MOSFET  structure  is  particularly 
convenient  for  such  measurements. 

(iii)  Photoionization  Cross  Sections 

The  photoionization  cross  section  of  the  indium  center  has  been  determined 
by  measuring  the  time  constant  of  the  decay  under  illumination  between  the  two 
limiting  conductivity  states  of  the  MOSFET  structure.  These  results  are  shown 
in  Fig.  8  on  page  63  of  this  report.  In  making  some  of  these  measurements,  we 
have  employed  a  fairly  unique  arrangement  using  a  circularly  variable  filter 
and  globar  source. 

Such  measurements  could  also  be  readily  extended  to  other  impurity  centers, 
such  as  the  O.lleV  level  reported  by  Hughes  in  their  indium  doped  samples. 

Such  measurements  can  provide  an  important  optical  verification,  in  addition 
to  the  thermal  measurements,  of  these  type  of  levels. 


(iv)  Proposed  Technique  for  Determination  of  Carrier  Capture  Cross  Sections 

Capture  cross  sections  and  capture  coefficients  can  be  measured  using  the 
MOSFET  structure  by  employing  an  extension  of  techniques  used  previously  on 
p-n  junction  devices.  These  techniques  can  be  employed  to  determine  both  the 
"small  signal"  capture  coefficient  or  capture  cross  section  when  the  carrier 
is  nearly  in  thermal  equilibrium  with  the  lattice  and  the  capture  cross  section 
when  the  carriers  and  impurity  centers  are  under  the  influence  of  electric 
fields.  Both  situations  are  relevant  to  photoconductive  detectors  which 
operate  over  a  range  of  applied  voltage  and  electric  field  conditions.  In 
addition  relatively  high  field  conditions  may  exist  at  the  p+  contacts  formed 
on  the  photoconductor. 

Figure  2  shows  the  proposed  technique  for  the  determination  of  the  cap¬ 
ture  cross  sections  near  equilibrium.  The  impurity  centers  in  the  MOSFET 
surface  depletion  region  have  been  arranged  to  be  in  the  negative  charge  state 
and  the  initial  conductivity  of  the  MOSFET,  I^s  ,  is  low.  The  surface  is 
then  accumulated  for  short  periods  of  time  attracting  holes  to  the  surface. 

If  this  accumulation  time  period  is  short  then  only  a  fraction  of  the  negative 
impurity  centers  can  capture  holes; 

*acc  ~  ^accumulation  <<:  ^cp  '  P) 
and  on  each  accumulation  cycle 

d(NT)/dt  =  -cp  p  • 
and  since  dt  =  t,„„  then 

aCC 

d¥*T  ■  -cp  e  •  ‘acc 
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Technique  for  Determining  Hole  Capture  Cross 
Section  at  Negative  Acceptor  Center  using 
MOSFET  Device 


Knowing  the  hole  concentration  in  the  substrate  from  the  resistivity  of  the 
photoconductor  then  one  can  determine  the  capture  coefficient  cp  by  observing 
the  number  of  pulses  required  to  completely  fill  all  the  negative  centers 
with  holes  leaving  them  in  a  neutral  charge  state  and  consequently  the  MOSFET 
with  a  much  higher  conductivity. 

The  capture  of  both  minority  carriers,  electrons,  and  majority  carriers, 
holes,  can  be  observed  under  high  electric  field  conditions  by  a  direct 
application  of  techniques  used  previously  on  p-n  junctions,  except  here  cap¬ 
ture  will  be  observed  in  the  surface  depletion  or  space  charge  region  of  the 
MOSFET  rather  than  in  the  depletion  region  of  a  p-n  junction.  This  is  il¬ 
lustrated  in  Fig.  3.  The  impurity  centers,  either  indium  or  gallium,  will 
be  arranged  to  be  either  in  the  neutral  or  negative  charge  state.  In  the 
first  case  shown,  the  centers  are  in  the  neutral  charge  state  and  minority 
carriers  are  injected  from  the  substrate  by  illumination  with  light  which 
absorbed  in  the  substrate  generating  minority  carriers  or  electrons.  These 
electrons  diffuse  and  then  are  collected  and  swept  across  the  depletion 
region  and  can  be  captured  by  neutral  centers.  By  measuring  the  photocurrent, 
Iphoto’  an<*  observing  the  time  constant  associated  with  the  change  in  impurity 
charge  state  one  can  determine  the  capture  coefficient,  cn. 

In  the  second  case,  the  impurity  centers  have  been  arranged  to  be  in  the 
negative  charge  state  and  holes  injected  by  illumination  of  the  surface  with 
very  short  wavelength  light  which  is  strongly  absorbed.  Holes  then  drift 
across  the  depletion  region  and  can  be  captured  by  negative  centers,  the  change 
in  charge  state  of  the  impurity  centers  will  modulate  the  conductivity  of  the 
MOSFET.  By  again  observing  the  time  constant  associated  with  the  change  in 
conductivity  and  the  photocurrent,  Ip^oto»  one  can  determine  the  capture  co¬ 
efficient  cp  under  high  electric  field  conditions. 


Note  that  in  all  cases  discussed  the  capture  coefficient  determined 
DOES  NOT  REQUIRE  A  KNOWLEDGE  OF  THE  IMPURITY  CONCENTRATION.  In  this  respect 
the  above  techniques  are  distinctly  different  than  most  previously  applied 
techniques  which  also  required  a  knowledge  of  the  impurity  concentration  which 
at  best  was  known  only  on  the  basis  of  an  indirect  determination.  These  tech¬ 
niques  and  particularly  the  first  for  the  determination  of  the  hole  capture 
coefficient  under  near  equilibrium  conditions  should  enable  an  accurate  result 
to  be  obtained  without  the  uncertainties  of  other  less  direct  techniques. 


